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ABSTRACT 

Aims. To present MUFFIT, a new generic code optimized to retrieve the main stellar population parameters of galaxies in photometric 
multi-filter surveys, and check its reliability and feasibility with real galaxy data from the ALHAMBRA survey. 

Methods. Making use of an error-weighted -test, we compare the multi-filter fluxes of galaxies with the synthetic photometry of 
mixtures of two single stellar populations at different redshifts and extinctions, to provide the most likely range of stellar population 
parameters (mainly ages and metallicities), extinctions, redshifts, and stellar masses. To improve the diagnostic reliability, MUFFIT 
identifies and removes from the analysis those bands that are significantly affected by emission lines. The final parameters and their 
uncertainties are derived by a Monte Carlo method, using the individual photometric uncertainties in each band. Finally, we confront 
the accuracies, degeneracies, and reliability of MUFFIT using both simulated and real galaxies from ALHAMBRA, comparing with 
results from the literature. 

Results. MUFFIT is demonstrated to be a precise and reliable code to derive stellar population parameters of galaxies in ALHAM¬ 
BRA. Using as input the results from photometric-redshift codes, MUFFIT improves the photometric-redshift accuracy by ~ 10-20%. 
MUFFIT also detects nebular emissions in galaxies, providing physical information about their strengths. The stellar masses derived 
from MUFFIT show an excellent agreement with the COSMOS and SDSS values. In addition, the retrieved age-metallicity locus for 
a sample of z < 0.22 early-type galaxies in ALHAMBRA at different stellar mass bins are in very good agreement with the ones from 
SDSS spectroscopic diagnostics. Moreover, a one-to-one comparison between the redshifts, ages, metallicities, and stellar masses 
derived spectroscopically for SDSS and by MUFFIT for ALHAMBRA reveals good qualitative agreements in all the parameters, 
hence reinforcing the strengths of multi-filter galaxy data and optimized analysis techniques, like MUFFIT, to conduct reliable stellar 
population studies. 
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1. Introduction 


Studying the stellar content of galaxies is crucial to understand 
their star formation histories (SFH), what in turn provides us 
with valuable information on the possible evolutive paths since 
their formation at high redshift down to the present time. Despite 
the large efforts and advances achieved in this topic during the 
last decades, it still remains as one of the most challenging and 
promising ways to understand galaxy evolution. 

Early attempts to study the stellar content of early- 
type galaxies were ba sed on colo u rs, from wide and nar¬ 
row band photometry (iBau 

iMcClure & van den Berghlll968t lEato 1973I) . and on empiri¬ 
cal synthesis of the populations using as reference basis the ob¬ 
served colours of nearby early-types. These early methods can 
be considered as the pioneers of the current photo-spectral fitting 
techniques, which are the main topic of the present paper. The 
above methods were gradually displaced by techniques based in 
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more specific features (lEabeill 1 9'?^ iPritchedl 1 977h that were de¬ 
fined in narrow spectral ranges. 

The arrival of absorption line-strength indices t o study the 
stella r content of galaxies dBurstein et al.l Il984t lEaber et ^ 
Il98-5h brought a significant breakthrough in the field. In this 
front, it is worth noting the Lick system of indices (iGorgas et al.l 
ll993l:IWorthev et al.ll994bl) . which for the last decades has been 
the standard for most spe ctroscopic studies i n stellar popula- 


tions in the optical (e. g. Trager et al. 119981; 

J 0 r 2 ensen 19991: 

Kuntschner et al.l 20011 Thomas etakl 2005 

; Bernard! et al.l 

20061 ISanchez-Blazauez et al.ll2006a; iGorgas et al.ll2007h. The 


combination of a certain number of absorption lines mainly sen¬ 
sitive to age, e. g. the Balmer lines, or to the metallicity, as 
traced by certain elements such as Ee, Mg, Ti, C, Ca, Na, 
etc., were proven to be an efficient way to break, at least to 
some ext ent, the well know n degeneracy between these two pa¬ 
rameters (IWorthevI 1 1 994all . The way to measure these features 
is delicately chosen to be very sensitive to a parameter of in¬ 
terest, focusing its study in small spectral ranges. By construc- 
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tion, line-strength indices are quite insensitive to the influence 
of extinction, and by fine-tuning their definition or combining 
the sensitivities of different indices, some of them may end up 
bei ng almost independerit frotn other parameters, e. g. meta llic- 
ity (IVazdekis & Arimotol 19991: Cervantes &JVazdekisi2009l) and 
a-element overabundances (iThomas et alJl200^ . 

In the last fifteen years, the development of stellar libraries 
in spectral ranges other than the optical has driven the defini¬ 
tion of new indices that allowed to extend t his kind of studies to 
other regions with unexplored sensitivities (ICenarro et al.ll20d^ 
iMarmol-Oueralto et al.ll20()^ . In addition, the index system of 
reference in the optical spectral range has been revisited and im¬ 
proved (see e. g. IVazdekis et al.ll2M^ thanks to the availability 
of much better stellar libraries at much better spectral resolution. 

It was with the arrival o f improved ste llar libraries, such as 
CaT (ICenarro et ahl I2001all5), ELODI E ( Prugniel & SoubiranI 
1200 Ih. ST E LIB dLe Borgne et alJ ]2003h . INDO-US 


Ivaldes et al.l 2()()4h. iMartins et all (l2005h. 


_ _ _ _ _ a nd M ILES 

Sanchez-Blazauez et al.l l2006bl: ICenarro et alJ l2007h . and 
the consequent e volutionary stellar popu l ation synthe- 


2003 

Gonzalez Dels 

2009 

Vazde 

ds et al. 1; 

Vazdekis et al. 

12012h. t 


Bruzual & Chariot 20031 : ^zdekiset^ 




Maraston et al. 


full 

spectral energy distribution (SED) of galaxies appeared as an 
alternative to line-strength indices. SED-fitting can also be used 
to de ri ve several physical properties of galaxie s (|Mathis_eLaL 


2005 iKoleva et ahl l200l ICoelho et al.l 1200^ IWalcher et al. 
2011 ; iLiu et al.ll2013l) . In fact, there is a growing number of 


public codes specifically devoted to c arrying out SED-fittin g 
with differe nt procedures, e. g. hyperz ( Bolzonell a et al.ll20()0l) . 
Le PHARE (lAmouts et al.l2002l: lllbert et al.l2006li . STARLIGHT 


(ICid Eernandes et ^|2005li. STECKMAP (lOcvirketal 

VESPA dT oieiro et al.ll2007 ), ULySS ( Koleva et al.ll200' 


FAST 


dKriek et al.ll2009l). SEDfit (Isawickill2012h. 

Nowadays, there is an increasing numb er of present an d 
future mu lti-filter surveys, e. g . COMBO-17 ( Wolf_e^^l2003), 
MUSYC (IC.awiser et al.hoObh. COS MOS (IScoville et akiboOTli 


ALHA MBRA (iMoles etallBOO^. CLAS H (iPostman et ahl 
l2012h. SHARDS (iPerez-Gonzdlez et al.l 1201 3ll . J-PAS 
(iBenitez et al.l l2014t) . and J-PLUS (Cenarro et al. 2015, in 
prep.), each of them with a vast volume of high-quality multi¬ 
filter data. These kinds of surveys pursue diverse goals with a 
common feature: sampling the SEDs of galaxies using top-hat 
and/or broad-band filters that mainly cover the optical range. 
Owing to this configuration, the retrieved SEDs are half-way be¬ 
tween classical photometry and spectroscopy, being in practice 
like a low-resolution spectrum whose resolution depends on the 
filter system (e. g. R ~ 20 for ALHAMBRA; R ~ 50 for J-PAS). 
Although multi-filter observing techniques suffer from the 
lack of high spectral resolution, their advantages over standard 
spectroscopy are worth to be noted: (i) The galaxy samples of 
multi-filter surveys do not suffer from selection criteria other 
than the photometric depth in the detection band of the survey, 
as all the objects in the field of view are observed. Eor a fixed 
observational time and similar telescopes, this leads to much 
larger galaxy samples than in multi-object spectroscopy, where 
achieving multiplexities larger than ~ 1 000 is a challenge at 
present; (ii) Unlike standard spectroscopy, the SED of galaxies 
observed in multi-filter surveys does not suffer from the typical 
uncertainties in the flux calibration that lead to systematic 
colour terms, as the photometric calibration of each individual 
band is independent from the rest. This advantage is crucial, as 
it is the overall continuum of the stellar population what in most 


cases dominates the diagnostic with SED-fitting techniques; 
(iii) With similar telescopes, the depth of multi-filter surveys 
is usually much larger than in that of spectroscopic survey, as 
direct imaging is much more efficient that spectroscopy; and (iv) 
Multi-filter surveys provide spatially resolved photo-spectra, 
similar to an lEU technique, allowing us to perform 2D stellar 
population studies in galaxies whose apparent sizes are not 
dominated by the point spread function (PSE) of the system. 

It is therefore clear that multi-filter surveys open a profitable 
way to advance in our understanding of galaxy evolution by pro¬ 
viding complete and homogeneous sets of galaxy SEDs down 
to a certain magnitude depth. Although there are several SED- 
fitting codes available, to cope with the calibration particularities 
of multi-filter surveys (see e.g. Molino et al. 2014), and given the 
vast amount of high-quality photometric data already available 
in the literature, and still to come in the next years, in this paper 
we present MUFFIT (MUlti-Filter FITtingfor stellar population 
diagnostics), a code specifically designed for analysing the stel¬ 
lar content of galaxies with available multi-filter data. This paper 
is mainly aimed to describe the code and its functionalities, set 
the accuracy and typical uncertainties in the retrieved stellar pop¬ 
ulation parameters, and demonstrate its reliability confronting 
with already existing stellar population results in the literature. 
The development of MUFFIT has been performed within the 
framework of the ALHAMBRA survey (see Sect. [2]), so, even 
though the code is generic and can be easily employed for any 
kind of photometric system, many sections in this paper are par¬ 
ticularized for the ALHAMBRA dataset. This allows us to show 
the code performances on real galaxy data, which is ultimately 
the best sanity check for any stellar population code. Despite in 
this paper we use galaxy data from ALHAMBRA, it is not our 
intention to exploit scientifically the data set here. In the next 
papers of this series, we will provide and exploit the stellar pop¬ 
ulation parameters retrieved for the whole galaxy sample in the 
ALHAMBRA survey. 

This paper is organized as follows. Section presents a 
quick overview of the ALHAMBRA survey, i. e. the photomet¬ 
ric dataset employed to develop the present work. In Sect. [3l we 
summarize the main technical aspects carried out by our code, 
MUFFIT, as well as the processes to obtain photometric colour 
predictions from models of single stellar populations (SSP) and 
the Milky Way (MW) extinction corrections. We show the ac¬ 
curacy and reliability of the stellar population parameters re¬ 
trieved with our code, together with the uncertainties and de¬ 
generacies expected for ALHAMBRA data in Sect.H] Section|5] 
presents a comparison study of the results retrieved from AL¬ 
HAMBRA galaxy data using MUFFIT with previous studies, 
including spectroscopic ones, and data from the literature, hence 
testing the reliability of our results. Einally, we provide the sum¬ 
mary and conclusions of this research in Sect, [b] 

Throughout this paper we assume a ACDM cosmology with 
Hq — 1\ km s“\ fiM = 0.27, and Da = 0.73. 


2. The ALHAMBRA survey 

The stellar population code that we present in this paper is gener- 
ically designed for all types of multi-filter surveys. However, 
we make use of the data in the ALHAMBRA survejQ to prove 
and test the reliability of our techniques, as in fact this code 
will be employed to analyse the stellar population properties of 
ALHAMBRA galaxies in forthcoming papers (Diaz-Garcia et 
al. 2015; in prep.). Therefore, throughout this work,we mainly 

* http://www.alhambrasurvey.com 
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present results, from both simulations and real observations, that 
are based either on the ALHAMBRA data or on its technical 
setup. In the following paragraphs we present a short summary 
of the ALHAMBRA survey. 

The ALHAMBRA survey provides a photometric dataset of 
20 contiguous, medium-band (FWHM ~ 300 A), top-hat fil¬ 
ters, that cover the complete opti cal range AA 3 500-9 700 A 
(see lAparicio Villegas et al.l 1201 Ol for further details) over 8 
non-contiguous regions of the northern hemisphere, amount¬ 
ing a total area of 4 deg^ of the sky (including areas in com- 
mon with other co smological surveys such as COSMOS, see 
iMolino et al.ll2014l for further overlapping areas). All filters in 
the optical range have very steep side transmission slopes, close 
to zer o overlap in wavel ength, a flat top, and transmissions 80- 
95% (iMoles et al.ll2008h . The magnitude limit is otab ~ 24 (5- 
sigma, measured on 3") for the 14 filters ranging from 3 500 to 
7 700 A, decreasing smoothly in the six re ddest filters reachin g 
down to Otab ~ 21.5 in the reddest one (iMolino et al.ll2014li . 
which is centred at 9 550 A. The optical coverage is supple¬ 
mented with the standard NIR J, H, and filters which have 
a 50% detection efficiency depth (point-like sources, AB mag- 
nitudes) of 7 ~ 22.4, // ~ 21.3, and K^, ~ 20.0, analysed in 
ICristobal-Hornillos et al.l (1200^ . The ALHAMBRA filter se0 
is design ed to optimize the accuracy of photometric redshifts 
(photo-z, iBemtez et alJ l2009l) . but due to their characteristics, 
it also provides low-resolution photo-spectra composed of 23 
bands, corresponding to a resolving power R ~ 20 in the op¬ 
tical. All the observations were done under a quality criterion, 
seeing < 1.6" and airmass < 1.8, using the 3.5 m telescope in 
the Calar Alto ObservatorjO (CAHA) with two cameras, the im¬ 
ager LAICA in the optical range and Omega-2000 for the NIR 
filters. At present, all the ALHAMBRA fields have not been im¬ 
aged yet, being the effective area for this work 2.8 deg^ with a 
total on-target exposure time of ~ 700 h (~ 608 h were dedicated 
for the optical bands, and ~ 92 h for the NIR ones), although the 
rest of fields will be completed reaching the expected total area 
of 4 deg^. 

The ALHAMBRA Gold catalogue 0 (IMolino etakllml . 
hereafter Gold catalog ue, is the reference catalogue for this 
work. As explained in IMolino et alJ (l2014l) . synthetic ^8141^ 
images were created, as a linear combination of individual fil¬ 
ters, to be used for both detection and completeness purposes, 
emulating the 7^8141^ band of the Advanced Camera for Sur¬ 
veys (ACS) in the Hubble Space Telescope (HST). Therefore, 
the G old catalogue prov ides 23-1-1 photometric AB magni¬ 
tudes dOke & Gunnll 1983b and errors for ~ 95 000 bright galax¬ 
ies (17 < mf^i 4 w < 23), being complete up to = 23. 

Throughout this work, the synthetic 7^8141^ photometry is re¬ 
moved from the analysis. Due to the existence of a PSF vari¬ 
ability among different filters, the photometry is corrected of 
PSF and aperture effects. In addition, and for the specific AL¬ 
HAMBRA case, we add quadratically an extra uncertainty of 
~ 0.025 (AB magnitudes) in each photometric measurement to 
account for potential calibration issues or uncertainties. 

For fu rther details of the AL HAMBRA survey, we refer the 
readers to lMoles et alJ(l2008li and iMolino et al.l (l2014t) . 


^ http://svo2.cab.inta-csic.es/theory/fps3/ 

^ www.caha.es 

http: //cosmo .iaa.es/content/alhambra- gold-catalog 


3. The code 

Although there are many public codes devoted to carrying out 
SED fitting in many different ways, e. g. hyperz, STARLIGHT, 
ULySS, VESPA, Le PHARE, FAST, or SEDfit, we are perform¬ 
ing our own analysis techniques to retrieve stellar population 
parameters from photometric SEDs, specifically designed for 
analysing the stellar content of galaxies from the ALHAMBRA 
survey, but being generic and easily adaptable to any multi¬ 
photometric survey. Secondarily, there is an increasing num¬ 
ber of large-scale multi-filter surveys; e. g. ALHAMBRA, J- 
PLUS and J-PAS, SHARDS, CLASH, MUSYC, COSMOS, or 
COMBO-17; offering a huge amount of photometric data that 
we can exploit to study the evolution of galaxies, opening a new 
path to explore the stellar population of galaxies, overall at in¬ 
termediate and high redshifts. Although these photometric data 
are like low-resolution spectra, these techniques present remark¬ 
able advantages respect to spectroscopy: they can go deeper with 
a better flux calibration (the calibration of each filter is inde¬ 
pendent of the rest of them), we can study the stellar content in 
each resolution element (similar to lEU techniques) with one ex¬ 
posure, and we can work with larger galaxy samples; hence, it 
would be a pity not taking advantage of these studies and not 
exploiting all the opportunities that they offer. 

The collection of analysis techniques, routines, and other 
tools that we are performing, are collected under the code name 
MUFFIT {MUlti-Filter FITting in photometric surveys), which 
is written in Python language, and it is mainly focused on re¬ 
trieving the stellar populations of galaxies whose SEDs are dom¬ 
inated by their stellar content. 

This section is subdivided in two extended sections. On the 
one hand, we show in Sect. 13.11 the main ingredients or inputs 
required to develop the analysis. These preliminary elements are 
basically composed by the SSP models, the photometric system, 
and the selection of a dust extinction law to treat properly the 
impact of dust on the model SEDs. On the other hand, the per¬ 
forming of the code is described in Sect. 13.21 being emphasized 
the description of some specific tasks. 

In Eig.[T] we outline the main structure of the code by a brief 
flowchart that summarizes the main features of the followed pro¬ 
cesses to set constrains in the stellar populations. We caution 
that the purposes of the flowchart is to help the reader follow the 
development of Sects. 13.11 and 13.21 and support the schematic 
comprehension of the large number of stages. 

The reader who may be primarily interested either in the re¬ 
liability of the code or in the comparison with results retrieved 
from the literature, may skip this section to continue with the 
self-contained Sects. |4]and|5] 

3.1. Main ingredients of the steiiar popuiation code 

In this section we describe the main input ingredients and 
preparatory tasks that are considered for the development of the 
stellar population analysis code that is presented in this paper. In 
particular, our code requires an input set of reference SSP models 
(Sect. l3.Llb . the photometric system of the data to be analysed 
(Sect. inr^ . and a set of recipes to take into account the intrin¬ 
sic and Milky Way extinction (Sect. I3.L3I) . The redshifts of the 
target galaxy data can be managed as an input ingredient or an 
output of the code, as it is explained in Sect. l3.L4l The flowchart 
on the left side of Eig. [Uparticularly illustrates the main ingredi¬ 
ents and preliminary work carried out by the code before starting 
with the analysis of the data. 
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Fig. 1. Flowcharts of the photometric model predictions (left) and the analysis techniques (right). A more detailed explanation of each step can be 
found in Sects. [TL^I3.1.3l and l3.2l The dashed row marks where both processes are related. Flowchart symbols represent standard tasks: ovals, 
start/end of a process (red); arrows, the direction of logic flow in the process; parallelograms, input/ouput operation (cyan); diamonds, a decision 
or branch to be made (yellow); and rectangles, a processing step (green). 


3.1.1. The SSP models 

The code has been designed to use SSP models as input 
templates for the comparative analysis of the stellar popula- 
tions of galaxies. Curren tly the code is ready to account for 
iBruzual & Charlod (I2003L hereafter BC03fl and MIUSCAT0 
(IVazdekis et al. 1Toi2l; iRicciardelli et alJ 12012) SSP models, al¬ 
though any other SSP spectral dataset can be easily imple¬ 
mented. 

BC03 is perfectly suited for SED fitting given the large spec¬ 
tral coverage of the models, from 91 A to 160 pm, allowing 
us to cope with most kind of multi-filter galaxy data, irrespec¬ 
tive of the redshift. For the present work, we assume ages up 
to 14 Gyr and metallicities [Fe/H]= -1.65, -0.64, -0.33, 0.09, 
an d 0.55, Padova 1994 track s (for f urther details and references, 
see IBruzual & Charlotj|2003l) . and a lChabrie^ (l2003l) initial mass 
function (IMF). 

MIUSCAT provides a sample of SEDs with a spectral range 
AA 3 465-9469 A and a r esolution of FWHM ~ 2.5 A , almost 
constant with wavelength (iFalcon-Barroso et al.ll201 ih . Despite 
the great colour calibration of these models, its spectral range 
is not enough for galaxies at intermediate redshift and further, 
missing the observed AFHAMBRA colours in the UV range. 
For this purpose, we extend the lower end of MIUSCAT models 
up to 1 860 A (A. Vazdekis 2015, pri v. comm.), using the N ext 
Generation Spectral Fibrary (NGSF. fHeap & Findlen 200^ . In 
addition, we complement these models with their photometric 
predictions for J, H, and K, which are adapted to predict the AF¬ 
HAMBRA NIR bands. Throughout this work, we use the models 


^ http://bruzual.org/ 

® http://miles.iac.es/ 


up to 14.13 Gyr with metallicities [Fe/H]= -1.31,-0.71, -0.40, 
0.0, and 0.22. We assum e a Kroupa Universal-like initial mass 
function (lKrour)all2001h . despite of its universality b eing a cur¬ 
rent matter of debate (see, e. g. iFerreras et al.ll2013h . In future 
works, we will shed light on the systematic variation of the IMF 
for the more massive galaxies in the AFHAMBRA database. 

By construction, the code can also use not only any other 
set of SSP models, but also any other kind of reference tem¬ 
plate spectra, e.g. spectra of real galaxies, as long as their main 
stellar population parameters (age, metallicity, IMF, extinction, 
and over-abundances) are assigned by the user. Throughout this 
paper we do not present this possibility, but concentrate on the 
performance of the code on the basis of the two SSP model sets 
mentioned above. 

3.1.2. Photometric system and synthetic photometry 

For a proper comparison between input SSP models and galaxy 
data, it is essential to build a reliable estimation of the synthetic 
magnitudes (or integrated fluxes) of the SSP template models in 
the same photometric system of the galaxies that need to be anal¬ 
ysed. This is computed by convolving the SSP model or galaxy 
reference spectra with the response functions of the photomet¬ 
ric system. In addition to taking the empirical filter transmission 
curves into account, to obtain a reliable photometric prediction 
it is advisable to account for specific characteristics of the ob¬ 
serving conditions and the instrumental setup employed for the 
photometric observations of the galaxies to be analysed. Among 
others, e.g. the transmittance of the optical system and/or the 
sky absorption spectrum where the observations were taken. It 
is remarkable the wavelength dependence of the quantum effi- 


Article number, page 4 of[3l] 



















































































L. A. Diaz-Garcia et al.: MUFFIT : A MUlti-Filter FITting code for stellar population diagnostics 


ciency of charge-coupled devices (CCDs), being typically less 
sensitive in the bluer and redder ends. If not accounted for prop¬ 
erly, this effect modifies the effective wavelength of such filter 
bandpasses creating a fictitious colour term in the synthetic pho¬ 
tometry of the reference models. In Fig.|2] the response functions 
of the ALHAMBRA photometric system is presented. It consists 
of 20 optical bands (left-hand side) and the ALHAMBRA 7, H, 
and NIR-bands (right-hand side). In this figure, all the effects 
explained above are already embedded. 

We compute th e synthetic photome t ry foll owing the method¬ 
ology described i nJPickles & Denag ^ ( 20101). wh ich is based in 
the HST svnphotfl package and in lBesselll (l2005h . 

As current detectors are photon-counting detectors, the num¬ 
ber of photons detected across a pass-band X is 


< = J ^F,Rx(d)dd, (1) 

where is the spectrum to convolve and Rx(X) is the response 
function of the filter X (also called sensitivity function in some 
previous work). Normalizing Eq. [T] we get the weighted mean 
photon flux density. 


^ fAF,Rx(A)dA 
^ f ARx(A) dd 


( 2 ) 


Some catalogues provide photometry in AB magnitudes, de¬ 
fined as 


show that the main, broader spectral features are easily distin¬ 
guishable after convolving, emphasizing the power of the AL¬ 
HAMBRA photo-spectra as halfway between classical photom¬ 
etry and spectroscopy. 

For the ALHAMBRA specific case, and because of the con¬ 
figuration of LAICA, we compute four photometric databases 
for the optical bands, one per CCD, as there exists discrepan¬ 
cies among CCD sensitivities and each CCD has its own set of 
filters. For the NIR-filters J, H, and K^, we repeat this process 
taking the Omega-2000 configuration (only one detector plate). 
In both optical and NIR synthetic photometry, we take into ac¬ 
count the filter transmission curves, the quantum efficiency of 
every CCD/camera, the sky absorption spectrum at CAHA, and 
the reflectivity of the 3.5 m-telescope primary mirror with the 
transmittance of the optical system. 

Due to both the large number of input model parameters 
(ages, metallicities, extinctions, IMF slopes and redshifts) and 
the intermediate-high spectral resolution of current SSP models, 
in general it is more efficient to build up our set of convolved 
models once at the beginning, rather than recomputing the model 
synthetic photometry every time the code is run. After comput¬ 
ing the synthetic photometry of all models, the photometric pre¬ 
dictions (fluxes and magnitudes) along with the main character¬ 
istics of each model are stored in a structured query language 
(SQL) database. 

A straightforward flowchart of the process to estimate pho¬ 
tometric predictions is shown on the left hand-side of Fig. [T] 


'«AB 


2-5 logio/v -48.6 , 


(3) 


where fy is the flux in ergs cm“^ Hz ' s '. To transform the 
weighted mean photon flux density into AB magnitudes, we 
compute the magnitude of the flux in the S TMAG system (sys¬ 
tem for calibrating FIST stars, IStonel[l99^ . which can be easily 
transformed to the AB magnitude system (Eq.|5]l. This interme¬ 
diate step is necessary due to the fact that the weighted mean 
photon flux density is established per unit wavelength, whereas 
the AB magnitude system is given per unit frequency. The mag¬ 
nitude across the bandpass X in the STMAG system, rnsj^, and 
in the AB system, otab,x, are 


msT,x =-2.5 logioF’P''-21.1, 


'«AB,x - insT,x - 5 logjo Tpivot -H 18.692 , 


'ipivot “ 


f ARx(A) dA 
^ f Rx(A) dA/A 


3.1.3. Dust Extinction 

Stellar population diagnostic techniques based on SED fitting 
over a large spectral coverage, as in this case, require the red¬ 
dening by extinction to be thoroughly taken into account to avoid 
potential misinterpretations of the integrated colours of the pop¬ 
ulation, e. g. older ages or higher metallicities, as well as to de¬ 
rive reliable stellar masses. 


extinction curve (e. 

2 . Prevot etal.l 19841: Massal 19871: MathisI 

199nHCardelli etal.l 

19891 Calzetti et al.ll200nl; O’Donnelll 19941 

Fitznatrickl 11999), overall on the bluer parts where the dust 


(4) 


(5) 


where Tpivot is the source-independent pivot wavelength, which 
is defined as 


reddening is more complex. The dust extinction curve is 
well reproduced usi ng the parameter Ry = AvlE{B - V) 
(ICardelli et al.l 1 198^ . which varies between 2.2 and 5.8 de¬ 
pending on the environmental characteristics of the diffuse inter 
stellar medium (ISM). Although the values of Ry may change 
depending of the line-of-sight, throughout this work we as¬ 
sume that the value of this parameter is Ry - 3.1, which is 
the mean value in the diffuse ISM of the MW (ICardelli et al.l 


(6) 


tion laws in our code (iPrevot et al. 

19841 

Cardellietal.lll989t 

lFitzDatricklll999tlCalzetti et al.ll2000 

), throughout this work we 


To illustrate this. Fig. [3 shows an example of the synthetic 
photometry determined for a SSP from BC03 with the AL¬ 
HAMBRA filter set. The black line is the flux of a SSP at rest- 
frame with solar metallicity, intermediate-age (5 Gyr), Chabrier 
IMF, and without intrinsic extinction. The colour squares cor¬ 
respond to the spectrum synthetic photometry following the 
process explained above, centred at their effective wavelengths 
(deft = f ARx(A) dA/ J Rx(A) dd), and the horizontal bars rep¬ 
resent the FWHM of each filter. This example is also useful to 


choose the Fitzpatrick reddening law (iFitzoatrickl 1 1 9991) as it 
reproduces well the extinction observed for MW stars with a 
prefer red mean value around Ry - 3.1 (ISchlaflv & Finkbeined 
1201 Ih . 

For extragalactic objects, there are two main sources of ex¬ 
tinction to account for. On the one hand, the dust intrinsic to the 
observed galaxy, which is redshifted with the galaxy system. On 
the other hand, the observed SED is reddened by the foreground 
MW dust in the observer’s reference system. It is important to 
note that this local MW extinction cannot be corrected together 
with the intrinsic galaxy reddening as the emitted flux is red- 
shifted before being scattered by the dust in our galaxy. As we 


http: //WWW. stsci. edu/institute/software_hardware/stsdaspi2^phd»elow, we separately face both extinction effects. 
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Fig. 2. Response curves of the ALHAMBRA filter set for the CCD 1 in the optical range (LAICA camera; one colour from blue to red per band) 
together with the ALHAMBRA J, H, and filters (Omega-2000 camera; dark red) to make the model synthetic photometry. 



Fig. 3. Synthetic photometry of a SSP using the ALHAMBRA photometric system. The black line is the SSP flux, the colour squares are the 
expected pass-bands and the horizontal bars represent the FWHM of each filter. 


Following a given extinction law, the intrinsic extinction is 
applied to the SSP template models before they are red-shifted 
and convolved with the photometric system. Throughout this 
work the values of Ay range from 0.0 to 3.1 (in bins of 0.1 in 
the range 0.0-1.0, and in bins of 0.3 from 1.0-3.1). The intrinsic 
extinction can be added as 

F, - F,,o X lO-O-^^^ , (7) 

where F,io is the SSP-model/template flux at rest-frame, after 
adding extinction, and A,i is determined by the extinction law, 
which can be chosen by the user. Since it is not clear how Ry 
varies within a host galaxy and among different types of galaxies, 
we keep constant the value to Ry = 3.1, i. e. the mean value 
in the MW. This helps to avoid degeneracies and to reduce the 
number of free parameters, which is already very high and time 
consuming. In s pite of the differe nt reddening laws have intrinsic 
differences ( see lFitzpatrickI 1 9^ . we do not assume errors in the 
SSP template models owing to such uncertain ti es. 

We use the dust maps of ISchlegel et al.l (Il998h . hereafter 
SFD, in order to deal with the MW reddening in the line-of- 
sight of each galaxy in our sample. The SFD dust maps provide 


E{B - V) values in different positions of the sky by estimating 
the dust column density. These estimations were calibrated us¬ 
ing galaxies and assuming a standard reddening law, to infer the 
existence of galactic dust between the observer and the sources 
beyond the MW limits. As the spatial resolution of SFD is low, 
FWHM ~ 6.1' and pixel size (2.372')^, MUFFIT makes a bilin¬ 
ear interpolation in the E{B - V) grid for every position of the 
target galaxies. 

MUFFIT applies a foreground extinction correction for each 
individual galaxy photo-spectrum using an extinction law for a 
value of E{B - V) and Ry. The most simple way to de-redden the 
photo-spectrum of a given galaxy is to compute the extinction in 
the effective wavelengths of the different Alters and then correct 
the source photometry using the equation 

Fa,c = F^,red X , (8) 

where Fa,c is the flux corrected of MW extinction for a given 
wavelength, Fai red is the observed flux (reddened), and A,i is the 
extinction factor given for a extinction law. Since the transmis¬ 
sion curves of the Alters are not completely flat and the shape 


Article number, page 6 ofl31l 









































L. A. Diaz-Garcia et al.: MUFFIT: A MUlti-Filter FITting code for stellar population diagnostics 


of the continuum is source dependent, this approximation may 
be inappropriate for those filters that exhibit a gradient in their 
transmission curves (e. g. the lower and upper end of the AL¬ 
HAMBRA optical bands, see Fig. la, especially in the spec¬ 
tral ranges where the observed spectrum is not flat. This effect 
would be interpre ted as a shift in the Alter effective wavelength 
(lFitzDatricklll999ll . and Anally, as a colour term in the spectral 
regions with strong gradients in flux, e. g. the 4000 A-break. To 
get a more reliable correction in this sense, the code carries out 
the de-reddening process of the data in three steps: 

- First, we pick up a set of models from BC03 (29 ages, from 
0.1 to 10 Gyr, four metallicities, [Fe/H] = -0.64, -0.33, 
0.09, and 0.56, and a Chabrier IMF) to be redshifted (redshift 
bin 0.01) and convolved with the survey photometric system. 
Before redshifting and computing the synthetic photometry, 
we add the intrinsic extinction {Ay from 0.0 to 1.0, in bins 
of 0.2 ) to the rest-frame BC03 models. Then, we carry out 
an error-weighted test to And the best fitting between the 
above models and the observed galaxy photometry. The aim 
of this step is not deriving physical parameters from the best 
fitting, but setting constraints on the shape of the continuum. 

- Secondly, we re-normalize the BC03 spectroscopic model 
associated to the best-fitting photo-spectrum. The synthetic 
photometry of this re-normalized model has to exactly re¬ 
produce all the observed photometric bands. 

- Finally, we apply Eq.0on the re-normalized model derived 
in the previous step, in order to obtain a de-reddened spec¬ 
trum that w e convolve with its related filter response-curve. 
We use the iFitzpatri^ d 19991) extinction laws to calculate 
A^, the value E{B - V) provided by SFD and Ry - 3.1, to 
de-redden all the galaxies of our sample. 

In particular, the iFitznatrickl (Il999h extinction law was built 
from the superposition of the extinction curves derived for a set 
of stars. Consequently, this extinction law contains intrinsic un¬ 
certainties, although we would accurately know the values of Ry 
and E{B -V). We account for the particular uncertainties of this 
law adding an error to the de-reddened photo metry of MW dust , 
F^ c, following the methodology explained in iFitzpatrickI (Il999h 
and assuming cru^ - 0. 

Cosmological fields, often the targets of multi-filter photo¬ 
metric surveys, use to be regions of the sky with low extinction 
values. In the particular case of ALHAMBRA, our main galaxy 
sample has MW extinction values of E{B - V) down to 0.04 
{Ay < 0.12 for Ry - 3.1) in all the cases. The colour term due to 
the MW dust in the ALHAMBRA survey may reach a maximum 
of Aotab ~ 0.15, and the stellar masses may be underestimated 
by 3% (8%) if we use the Ks (R) filter to estimate it. Although the 
stellar mass is not primarily affected by MW extinction in these 
fields, the colour term might change the retrieved stellar popula¬ 
tions and consequently the derived stellar mass (see Sect. l3.2.3l l. 
In ALHAMBRA there are no galaxies at low Galactic latitudes, 
I I < 5, where the MW temperature structu res are not duly 
resolved in the SFD maps (ISchlegel et al .11199^ . 

3.1.4. Redshifts 

The code is generically prepared to provide, together with the 
mass and the stellar population parameters of the galaxy, an es¬ 
timation of the photo-z. It is worth noting however that this code 
is not intended to be a photo-z code. Due to the large number 
of potential model parameters that the code plays with, when 
the redshift is set as a completely free parameter in the fitting 


process there exists a slight degeneracy with other parameters 
(like extinction; see Sect. I4.4b that tends to overestimate the de¬ 
rived photo-z. To overcome this effect, the code is also prepared 
to accept as initial constraint a list of redshift values for each 
target galaxy: either a list of nominal redshift values, hence the 
code only performs the fitting process at exactly these redshifts, 
or a complete probability distribution functions (PDF) of red¬ 
shifts, and then the code only accounts for the model redshifts 
within the PDF interval. Because of the good results we ob¬ 
tain, throughout this work we use as input redshift constraints 
the photo-z PDFs pr ovided by the ALH AMBRA Gold cata¬ 
logue using BPZ2.0 (iMolino et al.ll2014l) . It is noteworthy that 
the combination of our code with the ALHAMBRA photo-z con¬ 
straints further improves the quality of the input photo-z alone 
(see Sect. lO) . 


3.2. The core of the MUFFIT analysis techniques 

This section is devoted to the main technical features and pro¬ 
cesses carried out by our code in order to constrain the stel¬ 
lar population parameters of galaxies in multi-filter data sam¬ 
ples. We first describe in Sect. I3.2. li the way in which the 
minimization is performed, with the addition of a mixture of 
SSPs giving a remarkable improvement, specifically computed 
for each galaxy, in order to set more precise constrains in the stel¬ 
lar populations. In Sect. l3.2.2l we explain in detail the process to 
detect those bands that may be affected by strong emission lines, 
helping to understand the overall fitting process. Section [3. 2. 3I 
details how the stellar masses are calculated from the fittings. In 
addition, a Monte Carlo approach is performed to set constraints 
on the confidence intervals of the parameters provided by the 
code, detailed in Sect. I3.2.5I Finally, we describe how we man¬ 
age the k-corrections as result of the fittings in Sect. I3.2.6I The 
content of this section is outlined on the right-hand side of the 
flowchart (see Fig.[TJ. 


3.2.1. Ihex^ minimization and mixture of SSPs 

Our stellar population analysis technique is based on error- 
weighted ;t'^-tests between the multi-filter galaxy data and the 
template SSP models of choice. Since SSP models are usually 
normalized to a initial stellar mass and both the galaxy distance 
and its luminosity are uncertain in a general case (in fact these 
are parameters generally derived from the fit), it is required to 
add a normalization term, e, in the classical x^ equation. This 
term minimizes the value for every model-galaxy pair, be¬ 
ing the result only colour dependent. Our normalization way 
takes into account all observed bands and associated errors, be¬ 
ing more robust for multi-filter surveys as, at most, they only 
contain a few dozens of filters (e. g. 23 in ALHAMBRA). This 
way, all the meaningful filters contribute to determine the best 
solution of the fitting (without giving up one of the best bands 
in order to normalise), and there is no risk that the normalization 
band is affected by emission lines or cosmetic defects. 

Due to the fact that the number of reliable bands in each 
source may be different from one object to the other (for some 
objects, some filters may be rejected due to observational, cos¬ 
metic or calibration issues), in general we divide every by 
the number of available, safe filters in each case. Depending on 
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whether we are working with bandpass fluxes or magnitudes, the 
X~ definition can be expressed as 


Xn 


Nr ^ 


P X=1 L 


- (fim + mx) 


, and 


( 9 ) 



Wp 

Z 


- Sffx 


( 10 ) 


where Np is the number of available filters in an observed galaxy, 
(9™ * is the observed X-filter (magnitude or flux), cr™’^ its error, 
nix ifx) is the X-filter model prediction (single SSP or SSP mix¬ 
ture, more details below) and Sm (sf) is the normalization term. 
For our purposes, and Sf are written as 


r Np 
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y 

Zj f 2 
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( 12 ) 


which correspond respectively to minimizing Eqs. 0 and [TO] 
for each galaxy, i. e., - 0. As we will show later 

(Sect. 13.2.31) . by finding the best stellar population solutions for 
each galaxy, we can estimate its stellar mass from the e values. 

Note that Eq.|3(the equation used throughout this work) is 
assuming that the distribution of errors is Gaussian, when in gen¬ 
eral the distribution of magnitudes is not Gaussian since these 
are logarithmic measurements of flux. Erom certain signal-to- 
noise ratios, S/N^ 5 (or uncertainties cr“ ^ 0.22), the magnitude 
uncertainties are quasi-normally distributed, being this approach 
valid. Consequently, we encourage potential MUFFIT users to 
take fluxes instead of magnitudes when several galaxy bands are 
compromised by very low signal-to-noise ratios, S/N 4-5. It 
must be also taken into account that a certain minimum signal- 
to-noise ratio is required to determine reliable stellar population 
parameters without being dominated by degeneracies, as it will 
be shown later on in this paper. 

Once we have defined how to compute the fitting goodness, 
the next step is to compare our set of models to retrieve the most 
likely stellar population parameters. We carry out this process in 
two different steps. 


- First, we run the ^-test described above with the set of SSP 
models selected by the user (base models), making a first 
determination of the bands that may be affected by strong 
emission lines. In short, for every redshift step of the SSP 
models, the code looks for a flux excess in the galaxy SED 
with respect to the SSP model SEDs, for all those bands 
that could be affected by emission lines at the given redshift. 
A more extensive explanation on our technique of detection 
of emissions lines in multi-filter galaxy data is presented in 
Sect. 13.2.21 When this is the case, those bands potentially 
affected by emission lines are removed from the fitting pro¬ 
cess, and the ^-test is repeated again without the affected 
bands. In addition, rather than taking the parameters of the 
best SSP fitting, we carry out a Monte Carlo simulation us¬ 
ing the proper signal-to-noise ratios in each filter (further de¬ 
tails in Sect. l3.2.5l) . Prom the set of parameters retrieved dur¬ 
ing the Monte Carlo approach, we map the parameter space 
of compatible solutions (overall age, metallicity, extinction. 


o 

00 

00 

0) 

bfi 



0 2 4 6 8 10 

Agep, mock [Gyr] 


12 14 


Fig. 4. Empirical relation between the luminosity-weighted ages of 
mock galaxies made of random mixtures of two SSPs, AgeL, mock^ and 
the best age determination for such mock galaxies derived from a sin¬ 
gle SSP fitting, Agessp- The yellow curve is the AgeL, mock median for a 
given value of Agessp, and represents the typical offset in age that one 
may expect when interpreting the SED of a mixture of two SSPs by 
fitting a unique SSP. 


redshift, stellar mass and IMF), although at this stage we 
only focus on the retrieved distributions of age and redshift 
to carry out the next step: the mixture of two SSPs and its 
sub-sequent SED-fitting process. 

Secondly, according to the age and redshift distributions de¬ 
rived from the initial SSP analysis, we make a new database 
of models consisting of a mixture of two individual base 
SSP models. The mixture is computed only for the best red¬ 
shift solutions determined in the previous step. Eor each red¬ 
shift value, the two-model mixture is constrained to com¬ 
bine two SSPs, respectively younger and older than a certain 
age threshold, agex, that is related with the most likely age, 
agessp, inferred from the Monte Carlo analysis performed in 
the previous step. This is a reasonable assumption given that 
the stellar content of galaxies are usually th e result of com¬ 
plex SFHs with multiple s t ellar populations (Ferreras & Silld 
I 2 OOOI : iKavirai et alJl2007t iLonoce et al.ll2014ll . and the age 
solutions derived from comparisons with single SSPs can 
be considered, in first order, luminosity-weighted means of 
the ages of the individual, true populations. To determine the 
agex value that allows us to define the limit between younger 
and older SSP mixtures for each galaxy, we have studied 
the empirical relation between the luminosity-weighted ages 
of mock galaxies made of random mixtures of two SSPs, 
agCL, modi, and the best age determination for such mock 
galaxies derived from a single SSP fitting, agessp. In Fig- SI 
we present the result of this study. As expected, we ob¬ 
serve that agessp underestimates the real age, in particular 
for ageL, mock ~ 6 Gyr. The yellow curve in Fig. 0 ] repre¬ 
sents agex as a function of agessp- Once the age threshold is 
established, we generate all the possible SSP combinations 
(younger and older than agex), including as a new degree 
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of freedom the stellar mass weight of each component. For a 
general case with n components per mixture, each magnitude 
in the band X of the new mixture model is expressed as 


21 .' 

21 . 


nix,mix 

= -2.5 log 

10 


(13) 
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/-th SSP model and a, is the relative flux-contribution of the 

n 

SSP model in the i-th component, with ^ a/ - 1 and 0 < 

1=1 

or,- < 1. Note that in our case, we are mixing two SSPs and 
consequently n = 2. 

After mixing the SSP models as explained above, the code 
searches again the best fitting solution, repeating the detec¬ 
tion of emission lines with the mixture of models as ex¬ 
plained in Sect. 13.2.^ As in the first step using a single SSP, 
we do not only provide the best solution but map the com¬ 
patible stellar-population parameters by a Monte Carlo ap¬ 
proach, treating properly the errors in each band. This pro¬ 
vides an extra advantage when carrying out a statistical treat¬ 
ment of the results. We devote Sect. l3.2.5] to explain in detail 
how we explore the compatible space of derived parameters 
for each galaxy. 

Notice that, with this method and two SSPs, one database of 
mixed SSPs is particularly created for each galaxy, being more 
adequate and realistic than a single SSP fitting. As shown above, 
for a non-parametric SFH this represents a substantial improve¬ 
ment with respect to using one SSP only, which is not able to 
reproduce the colour of an underlying main red population plus 
less massive and later events of star formation. The mixture of 
two populations is a reasonable compromise, to improve signif¬ 
icantly the reliability in the determination of the stellar pop¬ 
ulation parameters of mu l ti-filter galaxy data (iFerreras & SilkI 
l200Ct iKavirai et al.l 1200'^ Lxmoce^^ 2014t) . In fact, it has 
been demonstrated (e. g. Rogers et al.l201(ll) that the mixture of 
2 SSPs turns out to be the most reliable approach to describe 
the stellar populations of young early-type galaxies, as well as 
a very reasonable approach for older galaxies, in this latter case 
only slightly surpassed by the use of chemically enriched expo¬ 
nential models. So, the two SSP model fitting approach may be 
considered in general, as a reasonable method for analysing the 
stellar populations of most kind of galaxies in a consistent way. 
Moreover, given that MUFFIT does not impose constraints on 
the metallicities of the 2 SSP mixture, this can provide hints not 
only for age evolution but also for a metallicity build-up. That 
being said, future versions of MUFFIT will also account for the 
use of different sets of SSP or r-models for the best choice of the 
user. 


3.2.2. Emission lines 

Nebular emission lines appear frequently in the SEDs of galax¬ 
ies, even if these are dominated by the light contribution of 
their stellar content. In particular, dealing with multi-filter galaxy 
data, filters affected by emission lines may present a substantial 
excess in flux with respect to any combination of SSP models, 
as the later typically do not account for the nebular emission 
physics. To guarantee the accuracy and reliability of the stel¬ 
lar population parameters derived during the fitting process, it is 
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Fig. 5. Spectral fitting examples for galaxies from the ALHAMBRA 
survey using MUFFIT and the MIUSCAT SSP models. The galaxy 
photo-spectra and their errors are given in red, whereas the best-fitting 
models are given in blue. Top panel corresponds to a quiescent galaxy 
for which no emission lines are detected. Bottom panel illustrates the 
case for a star-forming galaxy for which MUFFIT detects three bands 
affected by emission lines, in yellow. The dashed black lines indicate 
the wavelengths for Hj8, H ff, and [S ii]. Photometric redshifts are given 
at the labels inside. 


crucial to detect and remove those bands that can be significantly 
affected by emission. Not only due to the fact that they are not 
comparable to SSP models, but also since filters contaminated 
by strong emission lines tend to exhibit much large luminosities, 
hence lower photonic errors, than the rest of bands, and dominate 
our error-weighted SED-fitting techniques (see Egs. InlandfTOli. 
On the other hand, it is worth reminding that the presence of 
strong emission lines may also provide a fruitful information, 
since they contribute to the restriction of the feasible redshift 
intervals of the galaxy. The redshift constraints due to nebular 
emissions are additionally considered during the analysis. 

The emission line detection process of our code is dependent 
on the specific photometric system of the galaxy sample, as it 
only accounts for those emission lines that are typically strong 
enough to affect the photometry of the given filter set. It is obvi¬ 
ous that the broader the spectral filter width, the larger the equiv¬ 
alent width (hereafter EW) of the line that may be potentially de¬ 
tected at a fixed signal-to-noise. In this sense, the code is initially 
fed with a list of target emission lines that depends particularly 
on the filter set, customizable by the user, with emission lines 
such as [On] TT3726,3729, [Out] TT4959,5007, H Balmer’s 
series, [Sii] /1/16717,6731, etc. Due to the design of MUFFIT, 
we can also provide a list of typical AGN emission lines to re¬ 
duce their effects in the fittings, but broad AGN lines may affect 
2 or more ALHAMBRA filters, hence the AGN emission line 
detection criteria in MUFFIT might fail or be inaccurate. It is 
very important to note that an excessive list of emission lines, 
mostly when they are spread in the large wavelength ranges, will 
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eventually derive wrong results, as some bands may be removed 
in excess. It is therefore advisable to restrict this list to the lines 
that can present a measurable excess in flux, which mainly de¬ 
pends on both the Alter widths and the line intensity. For this 
reason, some bands can be forced by the user to remain in the 
SED-fitting analysis irrespective of whether they can be poten¬ 
tially affected by emission lines. For instance, in ALHAMBRA 
we do not expect to be sensitive to any flux excess in the NIR 
bands due to the presence of emission lines, so they are never 
removed during the fitting process even if the code detects a flux 
excess in any of them. 

Once we specify in the code the emission line list, the emis¬ 
sion line detection process is carried out in two steps. First, tak¬ 
ing into account the model redshift, we fit our models (single 
SSP or SSP mixture) without all the bands that could be poten¬ 
tially affected by the specified emission lines, and explore the 
residuals of the best fitting. If the residuals of any of the poten¬ 
tially affected bands present an excess in flux/magnitude larger 
than a limit value provided by the code user, Aotel, and these 
residuals are deviated beyond a band-error factor, ctel, the bands 
are considered to be affected by emission lines and are removed 
from the fitting process. Both constraints Aotel and ctel are re¬ 
quired, the latter to assure that the excess in flux is not due to 
photometric uncertainties, and the first one to avoid removing 
those bands with tiny observational errors that present little dis¬ 
crepancies with respect to the models. Finally, we repeat the fit¬ 
ting without the bands identified in the previous step, getting a 
new set of reliable values, clean from nebular contributions. 

For the ALHAMBRA case, we use Ajmel = 0.1, as for lower 
contributions the affected bands do not affect significantly the 
SED-fitting results retrieved with MUFFIT. In addition, we set 
(Tel = 2.5 as a reasonable statistical threshold to detect emission 
lines over the noise. Figure |5] illustrates two SED-fitting exam¬ 
ples of two galaxies from ALHAMBRA. The top panel of Fig. |5] 
illustrates a fitting of a quiescent galaxy in ALHAMBRA with¬ 
out strong nebular emissions, whereas the bottom panel shows 
a galaxy for which MUFFIT detects that some bands may be 
affected by emission lines (yellow squares). The red curves rep¬ 
resent the observed photo-spectra, while the blue curve is the 
best-fitting model after the detection of emission lines process. 
The yellow squares are the bands where the influence of an emis¬ 
sion line is ticked, in this particular case H/I, Hq'-i-[Nii], and 
[S It]. The dashed black lines point out the wavelengths for H a, 
Hy6, and [S ii] at the galaxy photo-z. For this case, the detection 
of the emission lines particularly contribute to strongly constrain 
the redshift range. Despite the ALHAMBRA resolution (FWHM 
~ 300 A), we note that strong emission lines can modify the fit¬ 
ting results. In some cases, even HyS shows significant contribu¬ 
tions in the ALHAMBRA data set. 

Since we are providing both those bands that may be af¬ 
fected by strong emission lines and the residuals from the SED- 
fittings, we can easily estimate the flux excesses in order to pos¬ 
teriorly transform them to equivalent widths. The advantage of 
this method is that our best SED-flttings are mixtures of SSPs 
that already include the corresponding stellar absorptions, hence 
the residuals can be directly related to the absolute nebular emis¬ 
sion. The main limitation, in general, comes from the low resolu¬ 
tion of the data, as in many cases some Alters can be affected by 
more than one emission line, like e.g. Ha and [Nii]. Still, as it 
will be presented in Section 15721 this technique opens new paths 
for future work on emission-line galaxies with multi-filter data. 


3.2.3. Stellar masses 

As we explain in Section 13.2.11 both the normalization term e 
introduced in the minimization equation and the intrinsic lu¬ 
minosities of the two best fitted SSPs are directly related with 
the total stellar mass of the galaxy. SSP models are usually nor¬ 
malized to an initial stellar mass of IMq, but this decreases with 
time accounting for the evolution of the most massive stars. This 
effect is properly taken into account for determining the Anal 
galaxy mass by applying a correction term to each SSP, /cssp, 
which is usually provided by the models. 

Taking into account the above considerations, the total stellar 
mass, M*_t, of a mixture of n SSPs (for this work n - 2) can be 
expressed as 

n n 

M*,t = 'Yu ‘ ^ ^ X Yj ’ 

/=1 /=1 

n 

= gf X X ^ /c,; SSP o'i , (15) 

/=i 

where M* , is the stellar mass of each population in the mixture, 
fitn is the normalization term d efined in Eq . fTSl d]^ is the luminos¬ 
ity distance in cm units (see lHogdll999h . Ki, SSP is the relative 
stellar mass correction for the f-th component in the SSP mix¬ 
ture, and a, is the relative flux-contribution of the SSP model in 
the i-th component (see Eqs.[T3]and fT4li . Throughout this work, 
the derived stellar masses are quoted including stellar remnants 
through /f, SSP, but for a more general case this parameter may 
not include remnants. 

3.2.4. Stellar population parameters of the SSP mixture 

The stellar-population parameters of the mixture of SSPs are 
estimated from the parameters of each component in the mix¬ 
ture. This can be done in different ways, which mainly depend 
of the weights assigned to the parameters of the different com¬ 
ponents. The most common definitions, provided by MUFFIT 
and employed in this paper, are luminosity-weighted and mass- 
weighted. The latter provides a more realistic information since 
it accounts for the total mass of stars in each population, hence 
assigning larger weights to the more abundant or dominant stel¬ 
lar populations. However, these populations may have very dif¬ 
ferent luminosities. In this sense, luminosity-weighted param¬ 
eters are more representative of the populations that dominate 
the observed spectrum, since the galaxy SEDs are predominantly 
leaded by the brighter populations, even if they are not dominant 
in relative mass. 

Throughout this work, the luminosity-weighted and mass- 
weighted stellar population parameters of a mixture of n SSPs 
(for this work n - 2), pL and pu respectively, are defined from 
the stellar population parameters of each f-th component (p,; 
age, metallicity, extinction, IMF slope, or [a/Fe]) as 

n 

2 a,- X Li X Pi 


27 cti X Li 
i=\ 


Pm 




M* 




Pi = 


n 

27 a, X /c,, SSP X Pi 

i=l 

n 

27 on X /f;, SSP 
(=1 


(17) 
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where a, is the relative flux-contribution of the SSP model in the 
i-th component, /<•,, ssp is the relative stellar mass correction for 
the i-th component in the SSP mixture, and L, is the luminosity 
of the SSP model in the observed spectral range. Note that both 
definitions agree when the parameter value is the same in each 
component. 

3.2.5. Determining the space of best solutions 

Because of the well known degeneracies among stellar popula¬ 
tion parameters, it is essential to perform a reliable analysis of 
the possible solutions (as mixtures of two SSPs) for each galaxy 
according to the uncertainties of the data. For this reason, rather 
than providing only the best fitting solution for each galaxy (it 
is well known that the most likely parameters are not always the 
best-fitting model parameters), our code accounts for the photo¬ 
metric errors of the multi-filter galaxy data to provide a set of the 
best fitting solutions, hence providing a set of probable values of 
redshifts, stellar masses, extinctions, and stellar population pa¬ 
rameters (ages, metallicities and IMFs) for each object. These 
values can be ultimately averaged according to their weights and 
frequencies to derive the average final parameters assigned to 
each galaxy and their errors. In this section we explain the pro¬ 
cesses and applied criterion to carry out this analysis. 

The determination of the best solutions space is based on a 
Monte Carlo method that, using the proper signal-to-noise ratio 
of each filter, seeks to obtain which parameter values are compat¬ 
ible within the photometric errors of the data. Since photometric 
uncertainties usually follow a normal distribution (or Gaussian), 
we assume an independent Gaussian distribution in each filter, 
centred in the band flux/magnitude, and with a standard devia¬ 
tion equal to its photometric error. It is worth noting that each 
filter is observed and calibrated independently from the remain¬ 
ing ones, so the errors of different filters are not expected to be 
correlated. 

For each galaxy, on the basis of the above Gaussian error 
distributions for its multi-filter data, MUFFIT generates Monte 
Carlo simulations (the number of realizations is defined by the 
user), ending up with a set of multi-filter data realizations for 
the same galaxy, all of them compatible within the errors. Ide¬ 
ally, the next step would be to run the test individually for 
each realization of the galaxy using the complete set of mod¬ 
els, but this is extremely time-consuming as the code plays with 
million of models (for the present research: 21 ages, 5 metallici¬ 
ties, 18 extinctions, 1 IMFs, 300 redshifts, and solar [(r/Fe]) for 
each fitting. Instead, to speed up this computational process, for 
each galaxy we perform a preliminary selection of SSP and mix¬ 
ture models that can play an important role in the fitting given 
the specific SED and errors of the galaxy. This pre-selection of 
models is carried out as it follows: 


distributions of the real galaxy multi-filter data. The corre¬ 
sponding 10 000;tf^ values between these realizations and the 
BFM, namely xl^, represent the range ofx^ values that one 
would expect just due to the photometric uncertainties of the 
real galaxy data. Note that this range can be very different 
among different galaxies. In MUFFIT, the limiting plausible- 
value, ;t'phot’ value that encloses the 68.27% (a 

Gaussian Icr) of the cumulative distribution function of the 
Xm values. 

- iii) Finally, the sub-sample of possible solutions for a given 
galaxy SED is constituted by those ones that fulfil the crite- 
^ Xbfm +Tphof sub-sample is consequently re¬ 
stricted to those models whose colours are statistically com¬ 
patible within the galaxy photon-errors. 


This way, the set of compatible best solutions for each galaxy 
is determined by generating A^ni Monte Carlo realizations of the 
galaxy SED data (throughout this work A^ni = 100) according to 
their errors, and then running our;^'^ minimization test for each 
galaxy realization using as input the sub-sample of preselected 
models. In each realization, we get a new BEM whose param¬ 
eters are ultimately weighted (cu j) with its x^ value to provide 
the most likely stellar population parameters together with their 
errors. Eormally, 


CJi = 


2 1/TmCj 

J=l 


(18) 


Wm 

< P >= 2 OJjPj 

J=l 


(19) 


^<p> — 




Nm 

Y^ajj[<p>- 


j=i 



( 20 ) 


where < p > and cr<p> are, respectively, the average stellar pop¬ 
ulation parameters (age, metallicity, extinction, redshift, stellar 
mass, IME, and [a/Ee], in a general case) and their errors, and 
Pj are the stellar population parameters associated to each BEM 
in the Monte Carlo realization with ax^ value equal 

In addition, the essential stellar parameters of each BEM ob¬ 
tained in the Monte Carlo iterations are also provided. 

Einally, we remark that the uncertainties of the parameters 
retrieved in this stage comprise not only the main parameters of 
the models, like ages, metallicities and IMEs, but also the extinc¬ 
tion, the redshift (if it is the case, within the interval provided by 
an external photo-z code) and the stellar mass. 


- i) After having run our code for a certain galaxy SED, and 
having obtained the x^ values for all the possible mixture of 
two SSP models (x^i^), we take the x^i^ value of the best 
fitting model (hereafter BEM),;^^gpJ^, i. e. the mixture of two 
SSPs with the greatest probability to be the solution, which 
corresponds to the lowest value. 

- ii) Since the parameter space of best solutions depends not 
only on the filter photometric uncertainties but also on the 
shape of the SED, the next step is to determine, for each in¬ 
dividual galaxy SED, the range of plausible x^ values that 
are expected according to the set of photometric uncertain¬ 
ties. To do this, MUFFIT performs 10000 Monte Carlo re¬ 
alizations of the BEM bands according to the Gaussian error 


3.2.6. K-corrected luminosities 

Once we have computed the best fitting models, we end up with 
a combination of SSP models that reproduce the colours of the 
galaxy photometric SED. Hence, the luminosity of the galaxy, 
its absolute magnitudes at any band, and the mass-luminosity 
relation is estimated from exactly the same combination of SSP 
models taken at rest-frame. Note that, independently of the phys¬ 
ical parameters linked to the best combination of models, the 
k-correction is model-independent since it properly reproduces 
the colours of a galaxy SED at a given cosmological distance, 
as long as the redshift is well constrained. If we compute the 
magnitudes for the different bands following this method, the 
main parameter that determines the k-correction goodness of 
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fitting is the photo-z accuracy. Since the set of SSP models 
does not contain emission lines templates, and our code re¬ 
moves them automatically during the fitting process, the pro¬ 
vided ^-corrections/luminosities only contains rest-frame pre¬ 
dictions about the stellar continuum, not about the nebular con¬ 
tent. 

To determine rest-frame magnitudes with the corresponding 
errors, for each galaxy we take all the best-fitting models recov¬ 
ered in the Monte Carlo approach (see Sect. 13.2.51) . average them 
and provide the average rest-frame magnitudes/colours and their 
standard deviations, hence considering the uncertainties in the 
photometry thanks to the Monte Carlo approach. It is notewor¬ 
thy that, at low redshifts, the uncertainties of rest-frame mag¬ 
nitudes may be very high, as the apparent magnitudes depend 
of the luminosity distance (oc d^^), which diverges at z = 0. 
This suggests that, to have better k-corrections in the most local 
Universe, more accurate photo-z are necessary. Despite this, the 
colour terms among different filters are not so affected by this 
effect, as the major impact is in the source luminosity and not 
in the rest-frame colours. To minimize this effect, we provide 
a second k-correction, for which we study the variability of the 
colours respect to an anchor band. In short, once we have all the 
rest-frame models recovered during the Monte Carlo method, 
the anchor band is the one that presents the lowest variability at 
rest-frame. In ALHAMBRA, this anchor band is usually a band 
in the red optical part (higher signal-to-noise ratios). This ap¬ 
proach turns out to be very useful, e.g. in order to make reliable 
colour-magnitude diagrams (CMD) at low redshift. 


4. Intrinsic uncertainties and degeneracies with 
ALHAMBRA galaxy data 

After having presented the main technical aspects of our SED 
fitting code in Section [2 and before presenting a comparison 
study between our stellar population results and similar previous 
data from the literature (see Sect. 13), the goal of this section is 
to study the accuracy and reliability of the stellar population pa¬ 
rameters retrieved with our code. Since this strongly depends on 
the photometric system of the data under study, it is important to 
note that, along this section, all the tests and predictions about 
uncertainties, degeneracies, etc. are particularly performed for 
the ALHAMBRA filter system. 

Since the code presented in this paper is particularly suited 
for the study of the stellar populations of galaxies whose SEDs 
are dominated by their stellar content, we begin with Sect. 14.11 
to build the CMD of the ALHAMBRA galaxy data, allowing 
us to make a proper selection of our target of galaxies and to 
compare our results with those published in the literature (Sec¬ 
tion |3- In Sect. 14.21 we check how the intrinsic uncertainties in 
the photometry of the ALHAMBRA filter system affect the typi¬ 
cal errors of the derived parameters, using a set of mock galaxies 
with well known input parameters. Eurthermore, the impact that 
the uncertainties of the input ALHAMBRA photo-z. Sect. 14.31 
have on the derived stellar population parameters is analysed in 
Sect. 14.41 Einally, we quantify the expected degeneracies among 
the derived galactic parameters of typical red-sequence galaxies, 
for the ALHAMBRA photometric system and different signal- 
to-noise ratios. 

4. 1. Selection criteria of ALHAMBRA red sequence galaxies 

It is well known that the CMD of galaxies exhibits a bimodal 
distribution with two main populations, usually referred as the 


"red s equence" (hereafter RS) and the "blue cloud" fe ell et al.| 


l2^lBaldrv et al.ll2^ [paber et al.ll20^ lEritz et alj|20l4 . A 


reat fraction of RS galaxies is mainly co mposed of early-types 
Strateva et al.ll200i] : Tcassata et al.ll2007h . but since the RS def¬ 
inition is clearly based on the observed galaxy colours, there is 
also a fr action of star formin g dusty-galaxies that may lie on 
the RS (IWilliams et al.ll2()()^ . To break the degeneracies be¬ 
tween quenched galaxies and dusty star-forming galaxies, there 
exist colour-colour diagno stics using NIR bands (IWilliams et al.l 
I 2 OO 9 I : lArnouts et al.ll201^ . and even methods to split the CMD 
into three populations ("red ", "blue", and "gr een") by fitting to 
a set of SED type classes (lEritz et al.ll2014l) . Eor the aims of 
this work, we just follow the classical method of the CMD 
(iBell et al'l l2004t lEaber et aP I2007I : lEritz et aI1l2014ll . A more 
detailed study of the contamination of star-forming, reddened 
galaxies in the RS will be given in (Dfaz-Garcfa et al. 2015, in 
prep.). 

To build the sample of RS galaxies, we firstly choose 
all the galaxies from the Gold catalogue^ with a statistical 
STAR/GALAXY discriminator parameter lower or equal to 0.5 
(Stellar_flag < 0.5), and imaged with 70% photometric 
weight on the detection image (PercW > 0.7), to avoid photo¬ 
metric errors in the galaxies close to the image edges. Secondly, 
we apply our analysis techniques over the full sample of AL¬ 
HAMBRA galaxies, using the set of MIUSCAT SSP models and 
the photo-z predictions included in the Gold catalogue, to auto¬ 
matically get their k-corrections (see Sect. 13.2.61) . Erom the k- 
corrections and the stellar masses, we can easily estimate their 
absolute magnitudes, that together with the rest-frame colours 
compose the CMD. We note that our CMD does not change sig¬ 
nificantly if we use another set of models, e. g. BC03, instead of 
MIUSCAT. In fact, this method is roughly model-independent 
as we are reproducing the luminosity and colours of the galaxy 
through the best mixture of two SSP models, irrespective of their 
parameters, hence the key here is to have a well-constrained 
photo-z (see Sect. 14.41 ). 

The RS and the blue cloud appear clearly separated when 
the CMD is constructed usin g the Johnson-like filters U and 
V (I.Iohnson & Morg^ Il953l) . In our case, for simplicity we 
select the ALHAMBRA filters F365W and F582W, as these 
are the ones whose effective wavelengths are most similar to 
U and y, respectively. The CMD of the ALHAMBRA galax¬ 
ies based on the F365W and F582W filters is presented in 
Pig .|3 where redder and bluer colours indicate higher and lower 
galaxy densities respectively. Pollowing the equation provided in 
iBelLeLak (120 041), w hich is compatible with the relation obtained 
in Pritz et al.l (l2014l) . we define the RS as those galaxies redder 


than the following colour-magnitude relation: 


mF365w-mF5&2w - L15-0.3z-0.08(Mf582W“5 log/2-I-20), (21) 


where m and M indicate apparent and absolute magnitudes in the 
Vega system. By simply visual inspection, it is clear that Eq. |2T] 
illustrated in Pig |6] as a red dashed line, splits properly the RS 
from the blue cloud, which already constitutes a first order check 
about the goodness of the SED-fitting. 


4.2. Photon-noise uncertainties 

To analyse the intrinsic uncertainties in the derived stellar popu¬ 
lation parameters of the galaxies due to the photon-noise errors 
of the ALHAMBRA photometry, we create mock galaxies con¬ 
sisting of a mixture of two random SSPs, in which we add photon 

* http://cosmo.iaa.es/content/alhambra-gold-catalog 
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Fig. 6. CMD of the ALFIAMBRA galaxies at different redshift bins. The filters F365W and F582W, both in Vega magnitudes, are used as a proxy 
for the Johnson filters U and V. Redder and bluer colours indicate, respectively, regions of the CMD with higher and lower galaxy densities. The 
dashed red line is the colour-magnitude relation that splits the RS from the blue-cloud for the mean redshift of the bin, as expressed in Eg. 1211 


noise according to the sensitivity of the ALHAMBRA filters and 
to the SED of the mock galaxies. Note that, by construction, this 
test is rather representative of the performance of RS galaxies. 
After adding noise, we run our code in order to derive the stel¬ 
lar population parameters of these mock galaxies, treating them 
as observed galaxies, but for which we know the real values of 
their parameters. The comparison between the input and the out¬ 
put parameters, as a function of the signal-to-noise ratio of the 
filters, allows us to conclude on the main topic of this section. 

We take the extended version of the MIUSCAT models (see 
Sect. B.l.TI) with the Kroupa’s IMF to develop these simulations. 
After adding different extinction values (Fitzpatrick’s law, Ay 
values ranging from 0.0 to 0.8 in steps of 0.1) at different red- 
shifts (from 0 to 1 and a step of 0.01), we randomly mix two 
SSPs with a series of constraints: 

- The weight of the younger population shall not be larger than 
30% in mass, and its age shall not be larger than 4 Gyr. The 
mass limit is required to avoid t oo low luminosity-w eighted 
ages, unlikely in RS galaxies (iKavirai et al.l 1200^ . and to 
guarantee the presence of old galaxies at all redshifts. 

- The age of the random SSPs cannot be much older than the 
age of the Universe at that redshift. Since the ages of SSP 
models are a discrete set of values, we state the limit to the 
first model that surpasses the age of the universe at each red¬ 
shift. 

- The extinction of both SSPs is the same. Although this may 
not be necessarily the case in general, it is a reasonable as¬ 
sumption as we are studying integrated stellar populations, 
which translates to an average intrinsic extinction that affects 
the projected incoming light from different populations. 

To properly sample the galaxy mass range, we assign ran¬ 
dom stellar masses in the range 9.5 < logjQ(M*/Mo) < 12.5. 
We repeat this process 2 000 times per interval of redshift, from 
0 to 1 in bins of 0.2, getting 10000 mock galaxies. As we ex¬ 
plain below, we study the impact of the signal-to-noise ratio for 
three cases {S jN - 10, 20, and 50). In each case, we construct a 
new random sample of mock galaxies, being the total number of 
simulations 30 000. 

After having built the set of toy mock galaxies, it is important 
to model accurately the way in which the galaxies are seen by the 
ALHAMBRA photometric system. It is in this point where the 
ALHAMBRA configuration plays an important role. The AL¬ 
HAMBRA characteristics (see Fig. and Sect. 13.1.21) are such 
that the reddest bands are not so deep as the rest of the LAICA 
filters. On the other hand, the SED of typical RS galaxies, even 
the youngest ones, exhibit a clear flux drop in the blue region, 
with a prominent 4 000 A-break in the middle. Therefore, we 


cannot assume either that all the filters present the same signal- 
to-noise ratio or that the signal-to-noise ratio among filters does 
not depend on the redshift. To carry out a realistic simulation, we 
take all the galaxies from the ALHAMBRA Gold catalogue for 
which the best-fitting corresponds to a RS galaxy (see Sect. EU), 
and compute, for every galaxy, the signal-to-noise ratios in each 
filter relative to the F199W filter, which is on average the band 
with the maximum signal-to-noise ratio at any redshift. By re¬ 
peating this process in different redshift bins, we determine how 
the signal-to-noise ratio changes along the SED as a function of 
the signal-to-noise in the anchor band F199W. These curves are 
shown in Fig.|2l and they account for the effective throughput of 
the telescope plus camera system, and for the average SEDs of 
RS galaxies. Note that the signal-to-noise ratios of the reddest 
filters up to z = 1 are strongly affected by technical features of 
the survey (mainly the depth in these bands), whereas the bluest 
filters are also affected by the SED shape. From the curves in this 
figure, it is easy to see how the 4 000 A-break moves from blue 
to red wavelengths when the redshift increases. Interestingly, at 
larger redshifts, the signal-to-noise ratio of the bluest filters starts 
to grow indicating larger fluxes in these bands, probably due to 
the pr esence of young populations in the galaxy (iFerreras & SilkI 
l2000h . which are easily observable at larger redshifts. Regarding 
the NIR filters, we check that typical RS galaxies become redder, 
on average, when they are observed at larger redshifts. 

To study the impact of different signal-to-noise ratios on the 
derived stellar population parameters, we add noise to the mock 
galaxies, taking in each case the suitable signal-to-noise ratio 
curve depending of its redshift. We build three samples of 10 000 
mock galaxies, and in each sample we force that the mean signal- 
to-noise ratio per mock photo-spectrum is S IN - 10, 20, and 
50 respectively; that is, for a galaxy with S/N - 20 at red¬ 
shift 0.5, the mean signal-to-noise for the 23 filters is 20, but 
in the bluest filter S INfiesw ~ 6 and in the anchor band (max¬ 
imum) SINfjggw ~ 30. The values of S/N - 10, 20, and 50 
respectively correspond to median apparent magnitudes for the 
detection band of mfsi 4 w ~ 22.6, 21.4, and 19.8 (ALHAMBRA 
RS galaxies and AB magnitudes). For the anchor band F199W, 
these values are almost identical. In ALHAMBRA, typical er¬ 
rors in the zero points due to calibration issues are ~ 0.025 (AB 
magnitudes), that correspond to a signal-to-noise ratio of ~ 50. 
Furthermore, most (^ 80%) of our ALHAMBRA RS subsample 
has a mean signal-to-noise ratio larger than 10, whereby these 
values {SIN - 10, 20, and 50) are suitable for our simulations. 

Although for the mock galaxies we take models with 0.0 < 
Ay < 0.8 and 0.0 < z < 1.0, for the mock analysis we use 
SSP models with redshifts up to 1.2 and extinctions up to 1.0 to 
avoid border effects in the parameter estimation. Concerning the 
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Fig. 7. Typical signal-to-noise ratios per filter for real ALHAMBRA red sequence galaxies at different redshift bins. The signal-to-noise ratios are 
normalized to that of the F199W filter. The observed data points account for the effective throughput of the telescope plus camera system, and for 
the average SEDs of red sequence galaxies. 


age estimation, we use the same constraint as in the mocks, i. e. 
depending on the redshift, the oldest ages are not allowed. 

Figure[8]illustrates the comparison between the input param¬ 
eters of the mock galaxies and the output parameters retrieved 
with MUFFIT, for the case S/N - 20 and all redshifts. Left 
panels present one-to-one comparisons for the input and output 
photometric redshifts, extinctions, luminosity-weighted ages and 
metallicities, and stellar masses. Right panels illustrate the distri¬ 
butions of the differences between the input and output values in 
each case, fitted to a Gaussian function (in red) whose mean and 
RMS are therein indicated. In addition, Table[T]provides the typ¬ 
ical mean differences and their RMS for different redshift bins 
and S/N= 10, 20, and 50. 

As expected, overall there is a very good agreement between 
the input and output derived parameters. This is not surprising 
as we are analysing mock galaxies made of mixtures of two SSP 
models, with the same SSP models as input of our code. In this 
sense, this test must be considered as a lower limit to the parame¬ 
ters uncertainties that we can expect for the forthcoming analysis 
of ALHAMBRA galaxies, just due to the photon-noise photo¬ 
metric errors. As a matter of fact, the total errors in the derived 
parameters are expected to be larger, due to potential differences 
between the spectro-photometric systems of the ALHAMBRA 
data and the models, independently of the SSP models of choice. 
Note also that real galaxies may be affected by ISM emissions 
or AGNs, which modify their SEDs with respect to a classical 
mixture of SSPs. 

Looking at the stellar mass plot in Fig.[8j there seems to be a 
slight overestimation of the stellar mass. These cases correspond 
to galaxies with z $ 0.02, for which little variations in the red¬ 
shift cause big changes in the luminosity distance, and therefore, 
in the retrieved stellar mass (see Eq.fTSl). This result suggests that 
in the very local Universe, more accurate redshifts are required 
to provide reliable stellar masses using the analysis techniques 
explained above. Fortunately, the very few local galaxies in the 
ALHAMBRA survey have a very high signal-to-noise ratio as 
well, whereby this overestimation is negligible in our case. 

Another case that worth to be explained is the one of mock 
galaxies with low extinctions and low metallicities. According to 
Fig.0 we are getting on average larger values. However this is 


an artifact of the simulations since there are not lower values in 
our set of SSP models. The important result in these plots is that 
we are still retrieving the right trend in the parameters, despite 
the border effects in the parameter space. 

The results in Table[T]are divided in different redshift bins as 
old ages are not allowed at high redshifts. It is noteworthy that 
all the parameters are better determined at high-redshifts than at 
low-redshifts at the same mean signal-to-noise ratio. First, be¬ 
cause at higher redshifts the galaxy SEDs are sampled with an 
equivalent higher spectral resolution at rest frame, whereby both 
redshift and age, that are sensitive to the 4 000 A-break, are bet¬ 
ter established (and consequently, the rest of parameters as well). 
Also, at high redshift the range of possible ages is shorter, and 
in turn they are younger, with lower degeneracies that their older 
counterparts. Finally, the bluest part of the z ^ 0.5 SEDs have 
larger signal-to-noise ratios. These filters act as anchoring bands 
to constrain blue-sensitive parameters (like extinction or metal- 
licity). This growth in the signal may be due t o an underlying 
young and less massive population in the galaxy (iFerreras & SilkI 
I2OOOI) . that is not strong enough to contribute in the optical range, 
but that dominates the flux in the NUV rest-frame regime (being 
visible at z ^ 0.5 in ALHAMBRA), reinforcing the necessity of 
using two components in the fittings. 

To conclude, these simulations are key to give us an idea of 
the typical issues that may appear in this kind of studies and 
the uncertainties that we expect just due to photon-noise pho¬ 
tometric uncertainties. These results show that one can robustly 
explore the stellar populations of galaxies in the ALHAMBRA 
dataset by use of the MUFFIT code presented here. 


4.3. Photometric redshifts in the ALHAMBRA survey 

Although the main aim of our code is not to determine redshifts, 
it is very important to check whether, for a general case in which 
the galaxies do not have any redshift information, the code is 
self-sufficient to estimate photo-z properly, at least to some ex¬ 
tent. Otherwise, the derived galaxy parameters may be wrongly 
estimated. 
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Fig. 8. Comparison between the input parameters of mock ALFIAMBRA galaxies, consisting of random mixtures two SSP models, and the 
output parameters retrieved with MUFFIT. The mock galaxies (RS-like galaxies) have average signal-to-noise ratios per filter of 20, according 
to the typical signal-to-noise ratio distribution presented in Fig. [7] Left panels show, from top to bottom, the one-to-one comparisons in redshift, 
extinction, luminosity-weighted age, luminosity-weighted metallicity, and stellar mass. Redder and bluer colours indicate, respectively, regions 
with higher and lower solution densities. The red dashed line indicates the one-to-one relationship. Right panels present the distributions of the 
differences between the input and output values in each case, fitted to a Gaussian function (in red) whose mean and RMS are indicated within the 
box. 
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Table 1. Typical uncertainties in the determination of redshifts, extinctions, luminosity-weighted ages, luminosity-weighted metallicities, and 
stellar masses, expected from running our code on ALHAMBRA red sequence galaxies at different redshift bins and for different S/N (10, 20, 
and 50). The random errors in the parameters are given as the mean and RMS of the best Gaussian function that reproduces the distribution of the 
differences between the input and output parameter values, as illustrated in the right panels of Fig. [8] 


Parameters 

0.0 < z < 0.2 

0.2 < z < 0.4 

0.4 < z < 0.6 

0.6 < z < 0.8 

0.8 <z< 1.0 

SjN^ 10 

^phot 

Ay 

logio Agey [yr] 
[Fe/H]L 
logio Af* [Mq] 

0.01 + 0.03 
0.10 + 0.22 
-0.01+0.19 
-0.09 + 0.29 
0.10 + 0.28 

0.00 + 0.04 
0.08 + 0.20 
-0.03 + 0.17 
0.01+0.30 
0.02 + 0.12 

0.00 + 0.03 
0.04 + 0.17 
-0.03+0.14 
0.03 + 0.26 
0 .01+0.10 

0.00 + 0.02 
0.04 + 0.15 
0.00 + 0.11 
0.00 + 0.23 
0.01+0.08 

0.00 + 0.02 
0.02 + 0.14 
0.03 + 0.10 
0.03 + 0.22 
-0.01+0.06 

S/N = 20 

^phot 

Av 

logio Agey [yr] 
[Fe/H]L 

logio Af* [Mq] 

0.00 + 0.01 
0.06 + 0.15 
-0.03 + 0.14 
-0.06 + 0.20 
0.02 + 0.15 

0.00 + 0.01 
0.04 + 0.13 
-0.03 + 0.13 
-0.03 + 0.19 
0.01 + 0.09 

0.00 + 0.01 
0 .01+0.10 
-0.01 +0.11 
-0.01+0.16 
0.00 + 0.07 

0.00 + 0.01 
0.02 + 0.09 
0.00 + 0.09 
-0.01+0.14 
0.00 + 0.06 

0.00 + 0.01 
0.00 + 0.08 
0.01+0.07 
0.00 + 0.12 
-0.01+0.05 

S/N = 50 

^phot 

Av 

logio Agey [yr] 
[Fe/H]L 
logio A7* [Mq] 

0.00 + 0.00 
0.01+0.08 
-0.02 + 0.11 
-0.03 + 0.08 
0.00 + 0.08 

0.00 + 0.00 
0.00 + 0.06 
-0.02 + 0.09 
-0.02 + 0.08 
0.01 + 0.06 

0.00 + 0.00 
0.00 + 0.05 
0.00 + 0.08 
-0.02 + 0.07 
0.00 + 0.05 

0.00 + 0.00 
0.00 + 0.03 
0.00 + 0.06 
-0.02 + 0.06 
0.00 + 0.04 

0.00 + 0.00 
0.00 + 0.02 
0.00 + 0.05 
-0.01+0.05 
0.00 + 0.03 


To do this, we run our code on the sub-sample of RS 
galaxies with spect roscopic redshifts in ALHAMBRA, from 
iMolino et al.l (l2014h . to set the accuracy of the ALHAMBRA 
photo-z. This sub-sample is built by the publicly available data 
of the spectr oscopic surveys that overl ap with ALHAMB RA 
(zCOSMOS iLillv et aLlEM GROTH. TPavis et alJl2()()^ and 
GOODS-N. ICooper et alJl201 ih . amounting to ~ 900 RS galax¬ 
ies up to magnitude < 22.5. For the purpose of this test, 

we make use of the photo-z predictions provided by BPZ2.0 in 
the Gold catalogue. In addi tion, we also take the photo-z con¬ 
straints provided by EAZY (iBrammer et al.ll2()()^ . with the de¬ 
fault configurations and templates, to asses whether MUFFIT 
also works similarly when the input photometric redshifts come 
from an external source. 

To provide a numerical value about the quality and accuracy 
of the photo-z, we simultaneously use definitions, which can be 
more or less useful depending of our purposes, for b oth accu¬ 
racy and catastrophic outliers. IBrammer et al.l (l2008l) proposed 
the normalized median absolute deviation, ctnmad, as a measure¬ 
ment of the photo-z uncertainties, since it estimates the deviation 
of the photo-z distribution without being affected by catastrophic 
errors. It is defined as 


ctnmad = 1 -48 X median 


|Az — median{Az)\\ 


1 "b ^si 


■I- 


( 22 ) 


where Az = Zphot - ^spec- Furthermore, we provide the RMS of 
the distribution Az/(1 +Zspec) that, in the following, we denote as 
o'z/il +Zspec)- Additionally, we us e two definitions for the rate of 
photo-z catastrophic outliers, as in iMolino et al.l (I20l4) . formally 
expressed as 


Ui 


m 


Az 


1 "b Zspec 

Az 


> 0.2, and 


1 +Z; 


> 5 X ctnmad- 


(23) 


(24) 


spec 
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Table 2. Quality of the photo-z retrieved for a subsample of RS galaxies 
in ALHAMBRA, when different method are applied: our code alone, 
BPZ2.0, EAZY, and our code using BPZ2.0 and EAZY as input values 
of redshift. Details on the definition of the quality values are given in 
Eqs.[22 |23l andl24] 



CTnmAD 

0-z/(l +Zs) 

m 

m 

MUFFIT 

0.0157 

0.0105 

1.6% 

5.8% 

BPZ2.0 

0.0104 

0.0076 

0.9% 

7.7% 

FAZY 

0.0102 

0.0083 

0.8% 

4.0% 

MUFFIT + BPZ2.0 

0.0087 

0.0070 

0.9% 

6.3% 

MUFFIT + FAZY 

0.0092 

0.0071 

0.7% 

5.5% 


On the basis of the above equations. Table |2] presents the 
quality of the photo-z determined for the subsample of RS galax¬ 
ies under different methods. First, we analyse the reliability and 
accuracy of the photo-z derived from our own code alone, that 
is, not using any photo-z value as input to constrain the solu¬ 
tion. This case can be directly compared with the values di¬ 
rectly derived from the BPZ2.0 and FAZY codes, showing, as 
expected, that photo-z codes make a better job to determine red¬ 
shifts from scratch. In addition, we analyse in the same way 
the quality of the photo-z derived from our code when the red¬ 
shift PDFs of BPZ2.0 are used as input parameters. MUFFIT 
explores the plausible stellar population parameters, managing 
the photo-z as another free parameter inside the redshift range 
of choice for the user. For the present work, we use the 1 - cr 
PDF range provided in the Gold catalogue, hence being all the 
photo-z weights equal to 1 inside the provided range (1 - cr) and 
0 beyond this range. According to Eqs.|22] |23 and|24] we ob¬ 
tain (Tnmad ~ 0.0087, yielding a rate of catastrophic outliers 
rji - 0.97% and 772 = 6.26%. For this specific case, the resul¬ 
tant photo-z are compared with their spectroscopic counterparts 
in Fig. |9] Similarly, using the photo-z constrains from FAZY as 
an input to our code (also 1 - cr, for consistency in the compar- 
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ison), we find ctnmad = 0.0092 ,771 = 0.76%, and 772 = 5.51%. 
From these results we conclude that, at least for RS galaxies, our 
stellar population code improves the redshift accuracy of clas¬ 
sical photo-z codes, when these are used as input for our tech¬ 
nique. This can be explained as our method plays with mixtures 
of much larger numbers (million) of SSP models, allowing to 
perform flexible SED fittings, hence providing a fine-tuned, sec¬ 
ond order correction to the redshift values of photo-z codes. In 
addition, the number of catastrophic outliers, 771 and 772 , is also 
marginally decreased on average. Regarding the shifts between 
photo-z and spectroscopic values, these are statistically insignif¬ 
icant 0.002), as we see in Sect. 14.41 on our stellar population 
results. As we show in Sect. 14.41 devoted to setting constraints 
on the uncertainties of the stellar population parameters due to 
the photo-z uncertainties, it is more important for our stellar pop¬ 
ulation aims to minimize the number of outliers, rather than de¬ 
creasing a little (Tnmad- 


4 . 4 . Impact of the photometric-redshift uncertainties 

One of the most critical parameters to determine reliable stellar 
populations and masses for a galaxy is the redshift. If the redshift 
is unknown or uncertain, it is obvious that the rest of derived 
parameters will not be reliable either. In this section, we aim 
to quantify the impact of the typical redshift uncertainties on the 
reliability of the retrieved parameters, and the maximum redshift 
errors allowed to reach our goals. 

To answer these questions, we focus again on the ALHAM¬ 
BRA data. To determine the direct impact of redshift uncertain¬ 
ties over the rest of derived parameters, we compare the results 
obtained from our code for the spectroscopic sub-sample of RS 
galaxies (see Sect. 14.3b using the redshift PDFs of BPZ2.0 as an 
input, with the results that we obtain with our code for the same 
galaxy sample assuming exactly their spectroscopic redshifts. It 
is worth noting that, in contrast to the previous simulations (e. g. 
Sect. 14.2b . we are using real galaxies to estimate these uncertain¬ 
ties. 

Figure [10] summarizes the results of this test. Left panels 
present the one-to-one comparison for the redshift, extinction, 
luminosity-weighted age, luminosity-weighted metallicity, and 
stellar mass. To facilitate the visual interpretation of the com¬ 
parison, red dots indicate those galaxies for which there exist 
large discrepancies between the photometric and spectroscopic 
redshifts, hence symbolizing the catastrophic outliers ( 771 ), as de¬ 
fined in Sect. 14.31 whereby yellow dots are the outliers related to 
772 . This is clearly seen in the top-left panel of Fig. [T0| where 
the redshift values are compared. In all cases, to guide the eye, 
the red dashed line indicates the one-to-one relationship. The 
right panels show the histograms of the differences between the 
obtained results with and without the spectroscopic redshift con¬ 
straint. The solid red line represents the best fitting to a Gaus¬ 
sian function, with its mean and RMS indicated in the panels. 
From the data in Fig. (TOjand the RMS of the differences given 
in the right panels, we observe that, overall, the stellar popula¬ 
tion parameters present very minor changes due to typical red¬ 
shift uncertainties. On average, such uncertainties account for 
- 0.03, =0.18 Gyr (0.03 dex), - 0.04 dex, 

and cr^ = 0.03 dex, all negligible when we compare them 
with the uncertainties introduced by the typical photon-noise 
(see Sect. wm . As expected, catastrophic outliers (red and yel¬ 
low dots) exhibit a larger spread in most parameters. 

To conclude, the typical uncertainties on the redshifts present 
a negligible impact on the main stellar population parameters. 


except if the galaxy is a catastrophic outlier. The uncertainties 
on the ALHAMBRA photometry and the model systematics are 
more crucial uncertainties in the present work. 

4 . 5 . Degeneracies 

In addition to estimating the uncertainties, it is crucial to know 
which kind of degeneracies may alter our results, to avoid un¬ 
veiling a finding that really is a degeneracy aftermath. To inter¬ 
pret the output properly, we must keep these degeneracies under 
control, reducing their impact as much as possible. 

Unlike stellar population diagnostic techniques based on 
local absorption features, e. g. classical line strength indices, 
whi ch also present the well known age-metallicity degener¬ 
acy dWorthevlllQ^ . our multi-filter stellar population code is 
colour dependent, as it tries to reproduce the galaxy SED by 
mixing SSPs over a wide wavelength range. Therefore, consid¬ 
ering only the age-metallicity degeneracy may not be enough, as 
we must evaluate any parameter that can modify the colour in a 
wide wavelength-range, i.e. we also have to include the intrinsic 
extinction as another degenerated parameter in our analysis. As 
we mention above, in this paper we assume an universal IMF. 
Otherwise, this parameter should be considered for the degen¬ 
eracies as well, as bottom-heavy IMFs exhibit redder colors than 
top-heavy ones. Moreover, since the degeneracies among param¬ 
eters strongly depend on the number and width of the Alters, the 
total spectral coverage, etc., it is worth noting that the results 
presented in this section only apply to the use of our code on 
ALHAMBRA data. 

To address the degeneracy issues in the most realistic way, 
we take as targets all the ALHAMBRA RS galaxies (see 
Sect. 14.1b with mean signal-to-noise ratios of ~ 20 and in cer¬ 
tain ranges of age and metallicity. These results come up af¬ 
ter having run MUFFIT with the ALHAMBRA galaxies using 
the MIUSCAT SSP models. The degeneracy estimation is done 
by taking all the stellar population values recovered during the 
Monte Carlo approaches at different stellar population bins (de¬ 
tailed below), and then stacking each retrieved distribution to 
build a whole distribution per bin, getting distributions among 
pairs of parameters (age, metallicity, and extinction). We char¬ 
acterise each distribution by setting confidence ellipses (2D con¬ 
fidence intervals) that enclose the results provided during the 
Monte Carlo process. These ellipses are obtained by the covari¬ 
ance matrix of each distribution, and they allow to parametrise 
the degeneracies through two parameters: by the ellipticity, de¬ 
noted as e, and by 0, the angle between the X-axis and the ellipse 
semi-major axis. The angle 0 is determined by the eigenvectors 
of the covariance matrix, as well as the eigenvalues of the co- 
variance matrix determine the axis lengths. If e is close to zero, 
this implies that the degeneration between the two parameters is 
not very significant, irrespective of the value of 0. On the other 
hand, if 0 is a multiple of ir/l (lie on any of the two the axes), 
both parameters are uncorrelated, and consequently there is not 
any degeneracy between them. Consequently, both parameters 
are necessary to confirm whether there exists a degeneracy or 
not. In fact, we can quantify the level of degeneracy between pa¬ 
rameters via the Pearson’s correlation coefficient, which mathe¬ 
matically reflects both e and 0 effects. Formally, 


2 (u - x)(yi - y) 


f xy 


I /^2 


Y^ixi-xYY.iyi-y) 

1=1 /=i 


(25) 
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Fig. 9. Comparison of the photo-z retrieved with our code using the redshift PDFs of BPZ2.0 as input values. The data correspond to a sub-sample 
of RS galaxies from ALHAMBRA with spectroscopic redshifts in the literature. Left panel illustrates the one-to-one redshift comparison for 
every galaxy, with the dashed line being the one-to-one relation. The bottom-right panel presents the differences between the photo-z and their 
spectroscopic counterparts for each galaxy, normalized by (1+Zspcc). The top-right panel shows the obtained distribution Az/(1 + Zspcc), indicating 
the accuracy parameters and the rate of outliers (Eqs.|^ and l24b at the inner box. 


where Xi and yi denote the value pairs of the parameters (age, 
metallicity, extinction) with means x and y respectively. The 
closer is to 1 (to -1), the larger is the correlation (anti¬ 
correlation), i. e. the degeneracy between parameters; on the 
contrary, a value close to 0 (-0.1 ^ r„j 5 0.1) suggests that 
the parameters are uncorrelated, and there would not exist any 
degeneracy. 

Regarding the parameter ranges, we take bins in age of 0.5 < 
AgeL < 1.0, 3.0 < AgeL < 4.5, and 7.0 < AgeL < 10.0 Gyr, 
whereas for metallicity we take -0.8 < [Fe/H]L < -0.6, 
-0.4 < [Fe/H]L < -0.2, and -0.1 < [Fe/H]L <0.1. These bins 
have been chosen to evaluate how the degeneracies vary along 
different stellar population parameters. Note that we do not es¬ 
tablish any extinction bin, because we have previously checked 
that the degeneracies at different extinction bins present negligi¬ 
ble differences in e and 0. We have also studied whether different 
redshifts can alter the degeneracy effects, as the redshift deter¬ 
mines the observed spectral range of the ALHAMBRA SEDs. 
We find that at higher redshift some degeneracies tend to de¬ 
crease, specially for young and low metallicity galaxies, but in 
general the degeneracies remains alike (less dispersion as well, 
see Sect. lO) . Thus, the three bins in age and in metallicity de¬ 
fine the nine intervals where the degeneracies of our parameters 
are explored. In each interval, we confront the age-metallicity, 
age-extinction and metallicity-extinction degeneracies. 

In Fig. [TT] we present the covariance error ellipses (in blue) 
that enclose, at the 95% confidence level, the distribution of the 
provided parameters during the Monte Carlo approach for all the 
ALHAMBRA RS galaxies in the nine age and metallicity ranges 
(see inner panels), and at redshift z < 0.4. The semi-minor and 
semi-major axes of each ellipse are shown in red, whereas their 
centres are represented with a yellow square. Table[2provides 0, 
e, and r^y for the same age and metallicity regimes. 


As intuitively expected, in all cases age and extinction are 
anti-correlated, in the sense that a reddening by extinction can 
mimic an older age, and viceversa. However, the behaviour of the 
age-metallicity and metallicity-extinction degeneracies is not so 
immediate. This is clearly a consequence of the role that extinc¬ 
tion plays in the analysis as a third freedom degree, absorbing 
partially the weight of metallicity in the classical age-metallicity 
degeneracy problem. Whilst older galaxies exhibit, as expected, 
clear anti-correlated age-metallicity and metallicity-extinction 
degeneracies, such anti-correlations may turn very mild or even 
positive correlations depending on the range of age and metal¬ 
licity. For instance, at the lower metallicity range, young and 
intermediate-age galaxies exhibit a positive degeneracy between 
age and metallicity, turning mild or negligible for young galax¬ 
ies with intermediate and high metallicities. Interestingly, young 
metal rich galaxies are essentially only subject to the extinction- 
age degeneracy. 

Finally, we have checked that the general degeneracy trends 
presented in this section do not vary qualitatively when com¬ 
puted on the basis of the BC03 models (see Table. O. For ages 
lower than ~ 1 Gyr, there is a clear degeneracy of both age and 
metallicity with extinction; and for ages older, there is also an 
age-metallicity degeneracy. 

5. Testing the perfomance of the code with 
ALHAMBRA gaiaxy data 

Once the technical details of MUFFIT have been described in 
detail, and the typical uncertainties and degeneracies among the 
derived parameters have been studied for the case of the AL¬ 
HAMBRA survey, in this section we apply the code to different 
sub-samples of galaxies in ALHAMBRA. The ultimate goal of 
this section is not to provide a thorough study of the stellar popu- 
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Fig. 10. Impact of a redshift uncertainty ~ 1% on the stellar population parameters. On the left, we present the comparative one-to-one of the 
obtained parameters without any constrain in the photo-z (Y-axis) versus the results forcing the redshift to its spectroscopic value (X-axis). 
The dashed red line represent the one-to-one relationship. Red, yellow, and blue dots are respectively galaxies for which Izspec - Zphotl > 0-2, 
0.044 < Izphot - ZspecI < 0.2, and Izphot - Zsped < 0.044. On the right, we have the histograms of the differences between the obtained results with and 
without the spectroscopic constraint. The solid red line is the best-fitting of the distribution to a Gaussian function, in the boxes there are showed 
both Gaussian mean and RMS. From top to bottom, we show redshift, extinction, age, metallicity, and stellar mass. 
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AlogipAgCL A[Fe/H]L AlogigAgeL A[Fe/H]L AlogipAgCL A[Fe/H]L 

- 0.5 0.0 0.5 - 0.5 0.0 0.5 - 0.5 0.0 0.5 - 0.5 0.0 0.5 - 0.5 0.0 0.5 - 0.5 0.0 0.5 



- 0.5 0.0 0.5 - 0.5 0.0 0.5 - 0.5 0.0 0.5 - 0.5 0.0 0.5 - 0.5 0.0 0.5 - 0.5 0.0 0.5 



- 0.5 0.0 0.5 - 0.5 0.0 0.5 - 0.5 0.0 0.5 - 0.5 0.0 0.5 - 0.5 0.0 0.5 - 0.5 0.0 0.5 



Fig. 11. Covariance error ellipses, at the 95% confidence level, of the stellar population parameters provided by the Monte Carlo approach in 
MUFFIT with respect to their most probable values, for different age and metallicity bins (see inner panels). Yellow squares indicate the ellipse 
centres and the red lines illustrate the minor and major axis of each ellipse. 


lations of ALHAMBRA galaxies, but rather to test the reliability 
of the stellar populations, emission line equivalent widths (here¬ 
after EW), stellar masses, extinctions and redshifts derived from 
our code in comparison with those published in previous work 
for either similar or identical galaxy samples. A forthcoming pa¬ 
per (Dfaz-Garcfa et al. 2015, in prep.) will present a complete 
analysis of the stellar populations of galaxies in ALHAMBRA 
making use of MUFFIT. 


5.1. Stellar masses and photo-z In the COSMOS survey 

Since the Fields in the ALHAMBRA survey partly overlaps 
with the COSMOS field, we can construct a sub-sample of RS 
galaxies (see Sect. I4. Il l in common between both surveys. After 
removing all the sources labe lled as stars in COSMOS (point¬ 
like sources, filbert et aDl2009l) and in ALHAMBRA, we end up 
with a sub-sample of 767 common galaxies up to redshift z < 
1 . 6 . 

This galaxy sub-sample has indeed an important added value 
for our testing goals, as it h as spectroscopic data in the zCOS- 
MOS lOk-bright catalogue (iLillv et al.l l2007l l2009l) . allowing 


Article number, page 20 of[3l] 
































































































L. A. Diaz-Garcia et al.: MUFFIT : A MUlti-Filter FITting code for stellar population diagnostics 


Table 3. Summary of the confidence error ellipses in Fig.[TT]for different parameter bins and MIUSCAT (top) and BC03 (bottom) models. 9 is 
the angle between the X-axis and the semi-major axis a in degrees (counter clockwise), whereas e is the ellipticity. To quantify the degeneracy 
between parameters, we provide the Pearson’s correlation coefficient, r^, where values close to 0 correspond to uncorrelated parameters. 
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for a calibration of the derived photo-z using the Le PFIARE 
code. In addition, using broad and medium b ands and fol l owing 
a SED-fitting technique with BC03 models, lllbert et ^ (1201 cil) 
estimated the stellar masses of the COSMOS galaxie s. They 
assumed a fixed Chabrier IMF (similar to the IMF of iKrouna 
120011) . a star formation rate ex e (with 0.1 < t < 30 Gyr), 
a unique sola r metal licity, an age grid of 0.1-14.5 Gyr, and the 
ICalzetti et al.l (l2000l) extinction law with 0.0 < E{B -V)< 0.5 
(large extinction values are only allowed for galaxies with high 
star formation). 

For the analysis of this sample we use MUFFIT, that imposes 
a mixture of two SSPs rather than an exponential star formation 
rate. We select MIUSCAT models with a Kroupa IMF (which 
slightly differs from the Chabrier IMF), for a wider range in 
metallicity ( -1.31 < [Fe/H] < 0.22). For our anal ysis we use the 
iFitzpatri^ (1 19991) extinction law rather than the ICalzetti et al.l 
(I 2 OOOI) law, as th e latter is oyer all more adequate for central star 
forming regions dCalzettill 19971) . 

In Fig.[T2]we present a one-to-one comparison of the photo-z 
and stellar masses obtained with our code on the AFHAMBRA 
data and those presented in the aboye work on the COSMOS data 
(lllbert et al.ll2009l 1201 Ol) for the sub-sample of 767 RS galaxies 
in common. Haying already constrained the reliability and accu¬ 
racy of the AFHAMBRA photo-z with respect to ayailable spec¬ 
troscopic redshifts (see Sect. 14.31 and 14.4b . in the top panels of 
Fig. [ 12 ] we confront our outcomes with the COSMOS photo-z, 


not only because of the qualitatiye similarity of both techniques, 
but also to check if there are systematics in our photo-z measure¬ 
ments that might cause any kind of systematic in the retrieyed 
AFHAMBRA stellar masses. From these plots we can see that 
both photo-z estimations are in yery good agreement, with one- 
to-one differences haying an RMS of just 0.015, and not finding 
statistically significant differences between both samples. The 
bottom panels of Fig.[T2]are deyoted to the stellar mass compar¬ 
ison. We find that the mean yalue of the one-to-one stellar mass 
differences is 0.04 dex, with a dispersion of 0.15 dex. 

As we explained in Sect l4.4l the typical uncertainties of our 
AFHAMBRA photo-z may haye an impact in the dispersion of 
the retrieyed stellar masses of up to 0.026 (see Fig. [loll, well be¬ 
low the obseryed yalue. On the other hand, the offset between 
masses cannot be completely explained by the mild offset be¬ 
tween redshifts (// = 0.005), as using Eq. [15] this difference im¬ 
plies a shift in mass % 0.015 dex. Preyious work haye already 
reported the non-negligible effect of using different model sets 
on the absolute yalues of the deriyed stellar masses (see e. g. 
iPozzetti et al.ll2007[ lllbert et ^l2010l) . To check this in our par¬ 
ticular case, we haye repeated the aboye analysis for the same 
sub-sample of galaxies with the BC03 models, instead of MIUS¬ 
CAT. The new mean difference between the stellar masses of 
AFHAMBRA and COSMOS galaxies is now -0.03 dex, with a 
similar dispersion of RMS = 0.15 dex. Consequently, the mild 
systematic between stellar masses is probably due to the SSP 
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Fig. 12. Differences between the photo-z (top panels) and stellar masses (bottom panels), computed with MUFFIT (Y-axis) and the values provided 
in COSMOS catalogues IX-axis Jllbert et alJ2009l . l201(]h . for a sub-sample of RS galaxies. On the right, we show the histograms of the differences, 
which are htted to a Gaussian distribution. 


model choice. Irrespective of the input set of SSP models, the 
dispersion between the stellar masses provided by COSMOS 
and the retrieved from ALHAMBRA data after running MUF¬ 
FIT remains RMS ~ 0.15 dex. This can be easily explained as 
the quadratic sum of the stellar mass uncertainties retrieved from 
both COSMOS and our proper technique in the common sample, 
being ~ 0.09 dex in COSMOS and ~ 0.11 dex in ALHAMBRA. 

As a general conclusion, despite the differences in the analy¬ 
sis techniques presented in this paper for ALHAMBRA galaxies 
and those performed by previous work for the COSMOS data, 
we find a remarkable agreement between the photo-z and the 
stellar masses derived for galaxies in common, showing the re¬ 
liability and robustness of our code. We confirm that the choice 
of SSP models and extinction laws has an impact in the absolute 
values of the derived stellar masses. In particular, we find that 
the stellar masses derived using the MIUSCAT models lead to 
mean values 0.07 dex larger that using BC03 models. 

5.2. Photometric EWs of emission iines 

During the SED fitting process, the detection and subsequent re¬ 
moval of the bands affected by nebular emission lines may be 
crucial to determine reliable properties of the underlying stel¬ 
lar content of the galaxy under study. The way in which the af¬ 
fected bands are detected and removed from the analysis is al¬ 
ready described in Sect l3.2.2l Here we analyse to which extent 
the emission line residuals retrieved from our fittings, based on 
photometric data, are reliable and still keep meaningful informa¬ 
tion on the true EWs of the nebular lines derived from classical 
spectroscopy. 

To build up a comparison galaxy sample, we first take all 
the ALHAMBRA galaxies that i) are in common with the 


data catalogues of the MPA/JHU 0 (hereafter MPA/JHU cata- 
logues) and ii) present nebular emission lines in their spectra 
(iBrinchmann et al.ll2004t iTremonti et alJl2004l) . This catalogue 
contains EWs and flux measurements of nebul ar lines for galax¬ 
ies in the SDSS DR7 (lAbazaiian et al.l l2009h . Such measure¬ 
ments already account for their corresponding underlying stel¬ 
lar absorptions, as they were calculated after the subtraction of 
appropriate SSP models. 

Since our stellar population code is focused towards the 
analysis of galaxies whose SEDs are dominated by their stellar 
content, for a fair comparison we systematically remove from 
the sample all the AGNs and QSOs (AGN, AGN_BROADLINE, 
QSO, and T2 types in SDSS), even if some of them could still be 
well interpreted by our code. In addition, galaxies with a signal- 
to-noise ratio lower than 5 in the EW continuum are removed, as 
well as those galaxies in the redshift ranges 0.112 < z < 0.114, 
0.123 < z < 0.125, and 0.146 < z < 0.148, to avoid EW contam¬ 
inations due to the sky line O i 45577, and all the galaxies larger 
than 4" (the SDSS fibers have 3" diameter) to minimize strong 
aperture effects in the photometry. Under the above constraints, 
there are 92 galaxies in common between ALHAMBRA and the 
MPA/JHU catalogues of SDSS. 

The detection and classification of emission lines in multi¬ 
filter surveys is clearly limited by the low spectral resolution of 
the data. Eor instance, at the ALHAMBRA resolution, and de¬ 
pending on the redshift, the emission line pairs H/l-[Oiii], Ha- 
[Nil], and even Ha-[Sii] can be unresolved because of their 
proximity in wavelength. To try to overcome this intrinsic lim¬ 
itation, rather than comparing the EWs of individual lines, we 
compare the total flux in excess along the observed spectral- 
range (AA 3500-9700 A for ALHAMBRA) with the total sum of 

^ www.mpa-garching.mpg.de/SDSS/ 
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the EWs measured in the MPA/JHU catalogues for the follow¬ 
ing strong nebular lines; H/3 /14861, [Oiii] A4959, [Oiii] 45007, 
Ho- 46563, [Nil] 46548, [Nii] 46584, [Sii] 46717, and [Sii] 
46731. We do not account for weaker lines as they might not 
be detected under the ALHAMBRA resolution. As explained 
in Section (3.2.21 we set the detection limits of emission lines 
in ALHAMBRA to a flux excess of Aiiiel = 0.1 and a signal- 
to-noise ratio with respect to the photometric error of the filter 
of (Tel - 2.5. Out of the 92 galaxies in common, there are 44 
galaxies in ALHAMBRA for which our code detects an emis¬ 
sion line in at least one filter, hence constituting the final galaxy 
sub-sample for the sake of LW comparisons. 

Since our SED-fitting technique is based on a minimiza¬ 
tion technique between the filter fluxes of models and real galax¬ 
ies over the full spectral range of the data, we must be aware of 
the fact that the best fitting solution may not match perfectly the 
local continuum around the emission lines, leading to random 
over/underestimations of the flux line. To minimize this effect 
for the sake of this study, given the Alter that contains the emis¬ 
sion line, fi, we define a local continuum for this band, as the 
mean value of the flux in the contiguous bands not affected by 
the emission line. Lor the model band that contains the modeled, 
corresponding stellar absorption, f/, we similarly define a contin¬ 
uum for the model, fj, in the same bands where is calculated. 
Lollowing the formalism for spectroscopic EWs, the equation for 
photometric EWs (in A) is 



where A4/ is the width of the band (~ 300 A for the optical 
Alters in ALHAMBRA) that contains the emission line, and the 
sum applies to all the bands affected by strong emission lines, X. 

In Lig.[T3]we present the photometric LWt, phot derived with 
our spectral-fitting techniques on ALHAMBRA data versus the 
spectroscopic LWt, spec computed from MPA/JHU catalogues, 
for the 44 galaxies in common that fulfilled the above selec¬ 
tion criteria. Blue dots indicate galaxies with only one strong 
emission line in their spectra, i. e. corresponding to a reason¬ 
able specific or individual emission line measurement in AL¬ 
HAMBRA, whereas red dots illustrate galaxies with more than 
one emission line. The black line is the one-to-one relationship. 
The dashed black line represents the minimum EWl, phot that we 
can detect in the ALHAMBRA data imposing Aotel = 0.1, for 
a unique emission line along the SED and a top hat Alter with 
LWHM~ 300 A. The red line is the linear regression to all the 
data points, considering both photometric and spectroscopic un¬ 
certainties. Overall, we obtain a good agreement between pho¬ 
tometric and spectroscopic EWs, where the bias between both 
measurements is // = 0.018 dex and RMS = 0.234 dex, demon¬ 
strating the feasibility and reliability of our code to determine 
the EWs of emission lines above a certain strength 13 A). In 
fact, our method can provide more robust determinations under 
certain conditions. This is the case, for instance, of the galaxy 
that deviates on the left side of the panel. The SDSS spectrum 
of this galaxy, at redshift z = 0.299, exhibits both Ha and H/J in 
emission. H/J emission is too weak to be detected in ALHAM¬ 
BRA with our criteria, whilst Ha, at 48526, falls in a very noisy 
region of the spectrum plenty of sky emission line residuals and 
telluric bands, which are hardly corrected in the SDSS data. This 
is not the case for the ALHAMBRA data, for which the contin¬ 
uum is better determined and where the absolute flux excess, for 
this particular case, becomes more reliable. 



Fig. 13. Comparison between the emission line EWs measured in galac¬ 
tic SDSS spectroscopic data taken from the MPA/JHU catalogues, 
FVFt, spec, and the emission line EWs derived from the residuals of our 
fitting code on the same galaxies in ALHAMBRA, EVLt, phot - Blue dots 
indicate galaxies with only one strong emission line in their spectra, i. e. 
with a reasonable specific or individual emission line measurement in 
ALHAMBRA, whereas red dots illustrate galaxies with more than one 
emission line. The black line is the one-to-one relationship, whereas the 
red line is a linear regression to all the points. The dashed hlack line is 
the minimum EVLt, phot that we can detect imposing Aiiiel = 0.1 for a 
unique emission line in one ALHAMBRA filter. See more details in the 
text. 


Linally we pay attention to the 48 galaxies for which our 
code does not detect any emission line in the ALHAMBRA data 
with the detection limits set in this paper. We find that 35 galax¬ 
ies (~ 75%) present logjoLWr, spec ^ 1-11 (~ 13 A), that cor¬ 
responds to our detection limit Aimel = 0.1, i. e. that indeed 
remain imperceptible under our detection constrain. About the 
remaining 13 galaxies, 10 of them have logjQLWr, spec > 1-11, 
but distributed along different lines, being all the lines individu¬ 
ally under the detection limit. Linally, for 3 galaxies we do not 
detect properly the emission lines for two reasons; first, since 
one of the emission lines (in this particular case Ha) is right in 
between two filters and the flux is split into both of them, not ful¬ 
filling the detection criteria in any of the filters, and secondly, for 
a wrong determination of the Zphot, that prevents the code from 
looking for emission lines in the right filters, besides the fact that 
a wrong redshift determination affects the quality of the derived 
continuum yielding a line residual under Aimel = 0.1. 

To conclude, we demonstrate that despite MUFFIT being 
mainly optimized for the analysis of the stellar populations of 
galaxies dominated by their stellar content, it is still reliable to 
detect and characterise the strength of strong emission lines un¬ 
der certain conditions that depend on the type of multi-filter data 
(e.g. filter width, signal-to-noise ratios per filter) we are work¬ 
ing with. In future versions of the code we expect to improve 
the algorithms of detection of emission lines with additional 
techniques and criteria (e. g, assuming intrinsic relations among 
lines), but this is out of the scope of the present paper. 
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Fig. 14. Spectral fitting of M32 as seen by ALHAMBRA using the 
MIUSCAT SSP SEDs as template models with the analysis explained 
above. The synthetic photometry of M32 is plotted in red, whereas the 
best fitting of a mixture of two SSPs to the spectrum of M32 is plotted 
in yellow. The bottom pannel shows the residuals of the best-fitting. 

5.3. The stellar populations of M32 

As a first step to test the reliability of the stellar populations de¬ 
rived with our code, we analyse the stellar content of M32, as this 
is one of the best known galaxies in terms of its resolved and un¬ 
resolved stellar populations. The spe ctrum of M32 used for this 
study is taken from the compilation of ISantos et al.l (l2002h . and it 
has been convolved with the ALHAMBRA filter set for the sake 
of this test, as if it had been observed in ALHAMBRA. Since 
the spectral range of this spectrum, AA 3500-10 000 A, is shorter 
than the filter coverage of ALHAMBRA, both the bluest optical 
filter and the three NIR filters have been rejected from the whole 
fitting procedure described in Sect. [3 In addition, to be able to 
explore the parameter space that is compatible with the best so¬ 
lution due to uncertainties in the photometry, and given that we 
have created a fake ALHAMBRA spectrum from a higher reso¬ 
lution spectrum, we add a synthetic error of cr ab = 0.025 in each 
filter, which is the expected error in the photometric calibration 
of ALHAMBRA (equivalent to a signal-to-noise ratio ~ 40). 

In Fig. [14] we present the complete spectrum of M32, in 
black, the M32 spectrum at the ALHAMBRA resolution, in red, 
and the best fitting derived from our code with a mixture of 
two MIUSCAT SSPs and a Kroupa IMF to the M32 ALHAM¬ 
BRA spectrum, in yellow. As consequence of the spectral range 
(AA 3 500-10000 A), the bluest and NIR ALHAMBRA filters 
have been rejected in the analysis. The obtained residuals are 
shown in the lower panel. 

It is clear from Fig. [T4| that the best fitting derived from 
our code reproduces well the observed spectrum at both low 
and high frequencies. The best fitting solution to a single SSP, 
as derived from our code in the first step, corresponds to a 
MIUSCAT model of 3.7 + 1.3 Gyr and around solar metal- 
licity ([Fe/H]= 0.02 + 0.14 dex). When the code is run com¬ 
pletely for the mixture of two SSPs, we obtain a luminosity- 
weighted age of 6.8 + 2.2 Gyr, a slightly sub-solar metallicity 
([Fe/H]=-0.08 + 0.14 dex), and extinction Ay = 0.28 + 0.08. 
Looking at the individual results for the two SSPs, we derive 
that the spectrum of M32 is well reproduced by an intermedi¬ 
ate age population of 2.1 + 0.5 Gyr and an older population of 


11.5 + 3.4 Gyr. We find that the weight on the stell ar mass of the 
young pop ulation is ~ 20% . Previ ous work, e. g. ICoelho et al.l 
(120091) and iMonachesi et aH (l2012h . obtain similar results, in the 
sense that M32 is not composed of a unique SSP of intermedi¬ 
ate age, but its stars were formed in at least two episodes of star 
formation, one ancient and the other one at intermediate ages. 

Note that, even though we are using 19 filters (instead of 23 
in a typical ALHAMBRA photo-spectrum), we still get a good 
agreement between the retrieved parameters and those derived 
by previous work making use of detailed spectroscopic studies, 
showing the great power of this kind of multi-filter surveys for 
stellar population studies. 

5 . 4 . Ages and metallicities of early-types in the local 
Universe 

Disentangling the stellar populations of early-type galaxies and 
their assembling histories is a key question for our understand¬ 
ing of galaxy evolution. However it is not our intention to ad¬ 
dress this point in this section. The aim of this section is only to 
explore the stellar content of a sub-sample of early-type galaxies 
in the nearby Universe from the ALHAMBRA survey, making 
use of MUFFIT, and compare our results with previous findings 
in the literature. Once again, this is an additional check to as¬ 
sess the reliability of the stellar populations derived from our 
techniques. In a forthcoming paper (Dfaz-Garcia et al. 2015, in 
prep.) we will carry out a more complete and systematic analysis 
of all the galaxies in ALHAMBRA, allowing to face these and 
other related questions. 

Our reference work is the paper by iGallazzi et al.l (120051 
in the following G05) and, in particular, the ages and metal¬ 
licities derived from spectroscopic analysis techniques for a 
sample of early-type galaxies located at z < 0.22 in SDSS. 
Their spectra were drawn from the S DSS DR4 (3" diameter 
fibres. [Adelman-McCarthv et al.ll20()^ spanning the full range 
of galaxy types (from actively star-forming to early-type galax¬ 
ies), covering the range AA 3 800-9 200 A with a resolution 
of R ~ 1 800, and Petrosian magnitudes in the r-band range 

14.5 < r < 17.77. To construct our sub-sample of early-types 
galaxies in ALHAMB RA we make use of the morphological 
catalogue provided by iPovic et ^ (l2013h . built using the code 
galSVM (desig ned to deal with low-resolution images at low and 
high redshifts, iHuertas-Comnanv et al.ll200^ that, following a 
Bayesian approach, classifies the galaxies morphologically. This 
catalogue contains more than ~ 1 500 early-type galaxies with 
redshifts down to z ^ 0.5 and a contamination lower than 10%, 
up to magnitude F613W < 22. To guarantee a fair comparison 
study, we only select from the above catalogue those early-type 
galaxies in ALHAMBRA with z < 0.22, i. e., in the same red- 
shift interval as in G05. With this constraint in redshift, we end 
up with a reliable sub-sample of ~ 400 early-type galaxies (mean 
signal-to-noise ratios S/N> 14 in all cases), in which a signifi¬ 
cant part (~ 65%) also reside in the RS. 

Following Figure 12 in G05, in FigllSjwe present the density 
contours of our results on luminosity-weighted ages, AgeL, and 
metallicities, [Fe/H]L, derived for our subsample of early-type 
galaxies in ALHAMBRA up to z = 0.22 using BC03 models 
and the photo-z constraints provided by the Gold catalogue. As 
in G05, the galaxy sample is split in stellar mass bins as indicated 
in the top labels of FiglTS] For each galaxy, rather than using 
the weighted values retrieved from the simulations, the whole 
set of results from the Monte Carlo simulations (see Sect. 13.2.5b 
are included in the plot. Darker colours correspond to the ages 
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Fig. 15. Luminosity-weighted ages and metallicities derived from MUFFIT, and the BC03 models, for a sub-sample of ALHAMBRA early-type 
galaxies at different stellar mass bins and up to z < 0.22. Purple crosses illustrate the semi-axes of the degeneracy ellipses for 1 - cr confidence 
level for the BC03 SSP models, as computed in Sect. 1431 in the same age and metallicites ranges. To guide the eye, yellow crosses illustrate an 
age of 6 Gyr with solar metallicity. 


and metallicity regions with higher population densities. Purple 
crosses illustrate the semi-axes of the degeneracy ellipses for 1- 
cr confidence level for the BC03 SSP models, as computed in 
Sect. 14.51 in the same age and metallicity ranges. 

From a purely comparative point of view, it is highly remark¬ 
able the similarity of the age and metallicity results derived from 
our code and the ones presented in Figure 12 of G05 on the basis 
of totally different spectroscopic analysis techniques. This is an 
excellent prove that the reliability of the stellar population tech¬ 
niques presented in this paper for analysing multi-filter galaxy 
data is comparable to the ones employed during the last decades 
for the analysis of spectroscopic data, like diagnostic diagrams 
based on line-strength indices, at least when they were applied 
to large volumes of data. 

From Fig[T5]we infer that low mass early-types, logjo M* 
10.3 dex, show a bimodal distribution in their stellar populations. 
There is a population of younger early-types with slightly lower 
metallicities that does not seem to exist at higher masses. On 
the other hand, there is a main population of older and, on aver¬ 
age, more metal-rich galaxies at the same stellar mass bin. The 
set of younger early-types of mild metallicities ma y be contam¬ 
inated by lenticular galaxies (iPoggianti et al.ll^Olh . and that are 
also found in G05. Using the rest-frame colours (see Sect. l3.2.6l l. 
we check that almost the totality of these "young" galaxies cor¬ 
responds to galaxies that reside in the blue-cloud, and proba¬ 
bly they are mainly composed of star-forming bulge-dominated 
galaxies. It is worth noticing that, despite this analysis is based 
on photometric data, the use of MUFFIT and ALHAMBRA al¬ 
lows us to be sensitive to this population and to characterise quite 
well their stellar populations. Not only obtaining a great agree¬ 
ment with spectroscopic studies, but also opening the possibil¬ 
ity of extending the mass limit up to lower stellar masses than 
typical spectroscopic surveys. We reinforce this result by repeat¬ 
ing the analysis using the MIUSCAT SSP models instead of the 
BC03, getting the same result. Since the detailed analysis of the 
stellar populations is out of the scope of this paper, this point will 
be addressed in a forthcoming work (Dfaz-Garcfa et al. 2015; in 
prep.). 

For intermediate stellar masses (10.3 log[gM*[Mo] ^ 
10.8), the "young" population tends to disappear, and conse¬ 


quently, the number of youngish and less metal-rich galaxies de¬ 
creases, being negligible for the higher stellar masses (10.8 $ 
logjo M*[Mo] $ 12.0), for which there is a clear predominance 
of old and metal-rich galaxies. For the lowest stellar masses, 
the spread in age and metallicity is apparently larger than for 
the most massive cases. Overall, our results suggest that mas¬ 
sive g alaxies are in average more metal-rich than less massive 
ones (iTremonti et al.l |2004 iGallazzi et alJ I2005L l2006h . hence 
the abundance of metals in a galaxy is related or linked with 
its stellar mass, showing a larger spread in the low mass end 
(also found in G05), except for the young metal-poor popula¬ 
tion. We also observe that the mean ages of massive early-types 
tends to be slightly older than their less massive counterparts, 
and, consequently, were formed at earlier epochs (higher red- 
shifts) than the lo w mass galaxies, in agreement with the "d own- 
sizing" scenario (ICowie et al.l 119961 : 1 Jimenez et al.l l2007h . The 
increase in the mean ages and metallicities for massive early- 
types was also found in G05 (equivalent results with very simi¬ 
lar age-metallicity relations using SDSS spectroscopy). Despite 
these stellar population differences are quite mild for the early- 
type galaxies (at least in comparison with late-type galaxies, see 
G05), it is worth remarking that MUFFIT, running on the AL¬ 
HAMBRA data, is still sensitive to the subtle changes in age 
and metallicity. On average, for the whole galaxy population in 
FigdSl we obtain that the increase in age from the low massive 
galaxies (9.0 log[QM*[Mo] ^ 10.3) up to the most massive 
ones (10.8 ^ logjg M*[Mo] ^ 12.0) is ~ 3 Gyr, being their mean 
ages 3 and 6 Gyr respectively. Similarly, the mean metallicity 
progressively increases from -0.35 to about solar metallicity. 

We perform several tests to assess whether the bi-modality in 
the populations of the less massive galaxies, 9.0 < logjgM* < 
10.3 dex, is driven by age-metallicity degeneracies. First, we ob¬ 
serve the degeneracies on the age-metallicity parameters for the 
whole sample of RS galaxies in ALHAMBRA at the parameter 
ranges of both distributions (Sect. 1431 see Fig.fTTIand Table |3ll. 
For both sub-populations (Age^ ~ 0.8 Gyr with [Fe/H]L ~ 
-0.7 dex; AgeL ~ 4 Gyr with solar metallicity), the degen¬ 
eracy contours do not present bi-modalities and they are well 
constrained by a unique ellipse (illustrated by purple crosses in 
Fig.fT3ll. In addition, we confirm that all the Monte Carlo realiza- 
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tions of each individual galaxy in the left panel of Fig.fTSlclearlv 
belongs to only one of the two galaxy sub-populations, evidenc¬ 
ing that degeneracies are not responsible for the younger and less 
metal rich population at the low mass regime. Moreover, we add 
noise (ctab = 0.05-0.20, corresponding to a signal-to-noise ratio 
S/N~ 20-5) to the high-mass galaxy sample and analyse them 
again with MUFFIT to see if there is any hint of bi-modality 
driven by degeneracies at low signal-to-noise regimes. Even in 
the worst case (S/N~ 5), less than 3% of the galaxies end up be¬ 
ing at the young and metal poor sub-population region, not ex¬ 
hibiting any bimodal pattern in the distribution. We conclude that 
there exists a true sub-population of "young" early-type galaxies 
in the stellar mass regime 9.0 < logjQ M* < 10.3 dex, not being 
a consequence of parameter degeneracies and the use of proba¬ 
bility distribution functions. 


5.5. Comparison with spectroscopic steiiar-popuiation 
studies 

A definitive step forward in the above analysis rests on the one- 
to-one comparison of spectroscopic galaxy ages and metallic- 
ities with the ones derived from MUFFIT. Interestingly, there 
is a sub-sample of galaxies in the MPA/JHU catalogues for 
which individual spectroscopic estimations of ages and metal- 
licities are provided (obtained following the methodology ex¬ 
plained in G05), and also imaged in the ALHAMBRA fields. 
G05 performed an age and metallicity diagnostic method based 
on a simultaneous fitting to 5 absorptio n line strength in¬ 
dices, most of them in the Lick system dGorgas et alJ Il99^ 
IWorthey et al. 199^ ), cons tituted by age-sensitive indices like 
D4000 (iBalogh et al.|[l9^ . HyS and H^A+HyA, and by m etal- 
sensitive indices like rMg,Le 1 (iBruzual & Charlol l2003h and 
[MgLe]’ (iThomas et al.l 120031) . the later weakly dependent on 
non-solar [a/Le] abundances. Each set of 5 spectral features is 
compared, through with the values provided by a set of 

models randomly generated from BC03, with different bursts of 
star formation and different fractions relative to the total stellar 
mass in a velocity dispersion range, to finally construct the PDL 
of the parameters being the weight of each model oc exp(-;t'“/2) 
(see details in G05). By crossmatching the ALHAMBRA galaxy 
catalogue with the above work, we find 80 RS galaxies (not spec¬ 
troscopically classified as either BROADLINE or AGN) in com¬ 
mon between both studies, with a mean signal-to-noise ratio per 
pixel in SDSS larger than 9 (under this constraint, the signal-to- 
noise ratio of the common ALHAMBRA galaxies is larger than 
18 in all cases). We establish this as a minimum threshold to 
obtain meaningful stellar population results from spectroscopic 
diagnostics based on line-strength indices. Despite G05 stated 
that metallicity is well constrained for spectra whose signal-to- 
noise ratio is larger than 20, this more permissive restriction in 
the SDSS signal-to-noise ratio increases the number of common 
galaxies, allowing us to explore the age compatibility of both 
methods for a larger number of RS galaxies, where the age ac¬ 
curacy is not so affected by the signal-to-noise ratio of the SDSS 
spectra (these details are extensively tested in G05). 

To keep the model consistency with G05, we feed MUFFIT 
with the SSP models of BC03 to explore, via ALHAMBRA data, 
the stellar content of these 80 RS galaxies in common. In addi¬ 
tion, to explore the impact of different SSP models in the re¬ 
trieved parameters, we repeat the same analysis with the MIUS- 
CAT SSP models instead (see at the end of this section). 

Pig. [I6] presents a one-to-one comparison of the spectro¬ 
scopic redshifts, luminosity-weighted ages, luminosity-weighted 
metallicities, and stellar masses given in the work by G05 for 


the sub-sample of 80 SDSS galaxies in common with ALHAM¬ 
BRA, and the photometric values determined from MUFFIT us¬ 
ing the BC03 SSP models for the same galaxies and the AL¬ 
HAMBRA data. In all panels, the dashed black line represents 
the one-to-one relationship, and the error bars at different pa¬ 
rameter ranges indicate the typical Icr uncertainties in the pa¬ 
rameters from both methods. Different colors indicate different 
extinction ranges (Ay < 0.3; 0.3 < Ay < 0.6; Ay > 0.6), as 
inferred from our code. In addition, to discuss aperture effects 
(SDSS spectra were taken in a 3" aperture, while ALHAMBRA 
photo-spectra are not restricted to a fixed aperture, which is de¬ 
termined by the synthetic band F814W), star-shape markers are 
assigned to galaxies with apertures below 4" in ALHAMBRA, 
which are, a priori, less affected by any potential aperture bias. 

Except for two galaxies, that have been removed from the 
plots since their ALHAMBRA photometry has been confirmed 
to be affected by nearby and bright stars, the spectroscopic red- 
shifts in Lig. [Tbb show an excellent agreement with our photo-z, 
with an RMS of ~ 0.008. As we expect from Sect. 14.31 we rule 
out that any difference in the stellar populations between the two 
sources can be due to uncertainties in the photo-z (see Sect. 14.4b . 

Concerning the luminosity-weighted age comparison, in 
Lig. fT6t > we find a good qualitative agreement between both 
methods given the uncertainties of both methods (see black 
crosses in Lig. [13, in the sense that lower (higher) spectroscopic 
ages correspond, respectively, with lower (higher) photometric 
ages from MUFFIT. Interestingly, according to MUFFIT, young 
galaxies also tend to be more reddened by dust than old galaxies. 
Lrom our results in this test-sample, the mean age (luminosity 
weighted) of the dusty galaxies (Ay > 0.6) is 1.8 Gyr, increasing 
to 3.9 Gyr for the intermediate extinction range, and rising up 
to 5.4 Gyr for galaxies for which we retrieve low dust contents 
(Ay < 0.3). This is an expected result, as it is well known that 
younger galaxies may still have remnants of gas and dust from 
recent star formation events, whereas older/quiescent galaxies 
use to have less dust content. We notice that this trend could 
be in principle explained by the age-extinction degeneracy (Sec¬ 
tion Ea). However, if this were the case we would not find any 
qualitative relation with the spectroscopic ages, as the line in¬ 
dices, by construction, are not significantly affected by extinc¬ 
tion. In addition, we emphasize that the galaxies in ALHAM¬ 
BRA for which MUFFIT retrieves young stellar populations, for 
ages down to 2.6 Gyr, are also classified as STARFORMING in 
SDSS. This also supports the idea that the retrieved extinctions 
with MUFFIT are very robust and they are not dominated by 
the degeneracy with age, as star-forming galaxies may present 
young populations with sig nificant amounts of dust. In fact, 
previous similar work (e. g. iLontana et al.ll20d5 iPozzetti et al.l 
120071: liibert et al.lbOlOh assumed in their codes that models with 
large extinctions are only allowed for star-forming galaxies. The 
key point out of this result is that MUFFIT is able to retrieve 
intrinsically the extinction of the stellar populations without as¬ 
suming any prior on the models or the galaxy type. Despite the 
qualitative agreement between MUFFIT and spectroscopic ages, 
with an RMS of ~ 1.6 Gyr, we also notice an offset between 
the ages of the two samples, with the spectroscopic ages being 
~ 2 Gyr older than the ones derived from our code for the AL¬ 
HAMBRA data. We discuss on possible reasons for this offset at 
the end of this section. 

In Lig.fTbb we can see that the metallicities present a quali¬ 
tative good agreement (despite the metallicity is more affected 
by uncertainties, see black crosses in Lig. [T6b . with an RMS 
of ~ 0.22 dex, although there is also a very small shift in the 
sense that our retrieved metallicities in ALHAMBRA tend to be 
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Fig. 16. Comparison of redshifts, ages, metallicities, and stellar masses between the spectroscopic study of iGallazzi et alJ i2005h . Y-axis, and the 
stellar populations retrieved from ALFIAMBRA with BC03 models, X-axis. In red, we present the galaxies for which our techniques stablish that 
may have large extinctions; in green, intermediate extinctions; and in blue, low extinctions. The dashed black line is the one-to-one relationship. 
The solid black line show the fitting between spectroscopic and photometric ages, accounting for the uncertainties in both measurements. The 
black crosses represent the average uncertainties of both techniques at different ranges. The star-shape markers are the galaxies with a radial size 
below 2". 


smaller (A[Fe/H] ~ -0.08 dex) than the spectroscopic ones de¬ 
rived for SDSS galaxies in G05. It is noticeable that the metal¬ 
licities of ALHAMBRA galaxies within an aperture of 4" (star- 
shape markers) present a better agreement with the spectroscopic 
measurements (A[Fe/H] 5 -0.05 dex) than the galaxies with 
larger apertures (dot markers; A[Fe/H] ~ -0.15 dex). As it 
is shown in G05, aperture effects and typical metallicity gradi¬ 
ents can lead to up to ~ 0.15-0.20 dex differences in metallic¬ 
ity for galaxies with ^ IO'^Mq. This reinforces the consistency 
and good agreement between both metallicity predictions (from 
SDSS and ALHAMBRA), since MUFFIT retrieves on average 
lower metallicities with respect to SDSS, and not larger, with a 
similar difference to those measured by G05 (~ 0.17 dex) ow¬ 
ing to possible aperture effects. Interestingly, as it is pointed out 
by G05, a signal-to-noise ratio larger than 20 is required for a 
reliable constraint of the metallicity in SDSS spectra. Unfortu¬ 
nately, the adopted signal-to-noise ratio limit substantially re¬ 
stricts our sample, excluding low luminosity galaxies with po¬ 
tential sub-solar metallicities. Moreover, in G05 it is also men¬ 
tioned that low-metallicity galaxies are more affected by uncer¬ 
tainties coming from their weak absorption features of [MgjFe] 
and [MgFe]’. Indeed, in our sub-sample, the lower the metallic¬ 
ity the larger the dispersion in the spectroscopic metallicity, as 
illustrated by the error bars. Instead, the metallicities provided 
by MUFFIT using ALHAMBRA data at the same regime are 
slightly better constrained, as, in the end, it is the overall stellar 
continuum what mainly determines the retrieved stellar popula¬ 
tions. 

Fig. [Tbl i exhibits a good agreement between the stellar 
masses of the two methods, with an RMS of ~ 0.19 dex, de¬ 
spite there is an offset of ~ 0.18 dex. This offset in the stellar 
mass can be mainly explained by the systematic differences of 
~ 2 Gyr between the ages of both methods, as this implies a shift 
in the mass-luminosity relation in the sense that older galaxies, 
at the same apparent magnitude, are also more massive galaxies. 

Before concluding, we aim at investigating on the potential 
origin of the ~ 2 Gyr offset in age derived in Fig. [Tbb . There 
are several potential reasons that could explain this offset: i) The 
age-extinction degeneracy; ii) The way to compute luminosity- 
weighted ages; iii) Aperture effects; and iv) Intrinsic systematic 
differences between both analysis techniques. 

i) Unlike SED fitting techniques, absorption line-strength 
indices are basically not sensitive to extinction, as they are 
defined in short wavelength ranges. If the ages derived from 
MUFFIT were severely affected by the age-extinction degener¬ 


acy, we would expect galaxies with very low extinction values, 
Ay < 0.05 to present a better agreement in the age compari¬ 
son of Fig. [T6t >. However this is not the case. By exploring the 
ages and metallicities of the galaxies with low extinction values 
{Ay < 0.05, according to MUFFIT and the ALHAMBRA data), 
checking that both metallicities, spectroscopic and photometric, 
remain in agreement without great differences, we obtain that 
the age difference is still ~ 2 Gyr. This discards the potential im¬ 
pact of the age-metallicity degeneracy, as well as the influence 
of using different extinction laws. 

ii) In G05, luminosity-weighted ages are computed accord¬ 
ing to the total r-band flux, whereas in MUFFIT this is done 
using the whole flux in all the ALHAMBRA bands. To explore 
whether this normalization difference could drive the age offset, 
we have recomputed the MUFFIT luminosity-weighted ages us¬ 
ing the ALHAMBRA F644W band, which has the most similar 
effective wavelength to the SDSS r-band. The results are essen¬ 
tially the same, hence not explaining the observed age offset. 

iii) It is well known that early-type galaxies may show radial 
variations of their stellar population properties, showing gradi¬ 
ents in metallicity and/or age. We already discussed above how 
the combination of different photometric apertures and the ex¬ 
istence of metallicity gradients has an impact on the metallicity 
comparison of Fig.fTbb). It is worth noting however that, in gen- 
eral, age gradients use to be shallower than metallicity gradients 

(IWu et al.l2005t[Sanchez-Blazauez et alj200^lLa Barbera et all 

l2012HEigenthaler & Zeilingetll2013ir so aperture effects are ex¬ 

pected to be smaller as well. Nevertheless, shallow gradients in 
all the parameters can also be found (e. g. Gonzalez-Delgado 
et al. 2015, in prep.). To assess this effect, we focus on galax¬ 
ies whose photometric apertures in ALHAMBRA are down to 
4" (star-shape markers in Eig. [T6l l. not far from the SDSS fiber 
aperture. We observe that the age offset is significant even for 
these galaxies, hence aperture effects are rejected to explain the 
age offset too. 

iv) After the negative results of the three previous tests, the 
existence of intrinsic systematic differences between the two 
methods seems to be the most plausible reason for the differ¬ 
ent absolute values of the derived ages. The discrepancies be¬ 
tween the analysis of spectral features versus colours, together 
with the assumptions of different SEHs (exponentially declin¬ 
ing tau models in G05, versus a mixture of young-Fold SSPs in 
this work), may be responsible for the age offset. To shed light 
in this last item, we aim at constraining a purely mathematical 
problem; the potential differences between the luminosity ages 
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Fig. 17. Comparison between the luminosity-weighted ages of a sub¬ 
set of RS T-models from SSAG and the ones derived by MUFFIT for 
the same models employing a mixture of 2 SSPs. The dashed black 
line represents the one-to-one relationship, while the red line is the sim¬ 
ple linear regression of the data points. The solid black line illustrates 
the fitting between the ages provided by G05 and the ones provided by 
MUFFIT using ALHAMBRA data (see Fig.IHr). 

derived from parametric r-models and the ones derived from a 
non-parametric mixture of two SSPs. To do this, we make use of 
the set of r-model s Synthetic Spectral Atlas of Galaxies (SSAG, 
iGladis et al.ll20l4 priv. comm.) very similar to the r-models em- 
pl oyed in G05. SSA G models are based in the recipes described 
in IChen et al.l (12012h and have been constructed following a ex¬ 
ponentially declining SFH and BC03 models, that may randomly 
suffer an instantaneous and random burst during different periods 
of time. The SSA G also includes intrinsi c extinctions, following 
the dust model of ICharlot & Falll (l2000l) . and different velocity 
dispersions. To create a sub-sample of RS galaxies, we select all 
the SSAG galaxy models whose colours satisfy U-V k. 2.0 (AB- 
system). After convolving SSAG models with the ALHAMBRA 
filter set, we run MUFFIT using as input the same BC03 models, 
in concordance with SSAG, and compare derived luminosity- 
weighted ages with the input SSAG ones. The result of the 
comparison is exhibited in Fig. [ni There appears a system¬ 
atic offset between ages of AAgeL ~ 1.8 + 1.7 Gyr (red line in 
Fig. [Ill), which fully explains (qualitatively and quantitatively) 
the ~ 2+1.6 Gyr offset (solid black line in Fig. fTTl) found in 
the previous comparison between the spectroscopic ages of G05 
and the ones retrieved using MUFFIT and ALHAMBRA data, 
as due to the mathematical differences in the diagnostic input 
models (mixture of SSP models vs exponentially declining SFH 
models). 

Finally, we have also investigated the impact of using a dif¬ 
ferent set of SSP models for the stellar population comparison. 
In Fig.[T8]we present the same comparison of Fig. [16] but in this 
case having used the MIUSCAT SSP models, instead of BC03, 
to analyse the ALHAMBRA data with MUFFIT. Except for a 
slightly larger difference in metallicity, in the sense that MIUS¬ 
CAT models tend to predict lower metallicities, the rest of pa¬ 


rameters compare similarly to Fig. [Tb] This is important to as¬ 
sess the impact of different SSP models in the absolute values of 
the derived parameters. 

6. Summary and conclusions 

The arrival of present and future large-scale multi-filter sur¬ 
veys (e.g. COMBO-17, COSMOS, ALHAMBRA, SHARDS, 
J-PLUS, and J-PAS) promises the availability of formidable 
datasets for many purposes in Cosmology and Astrophysics. 
These photometric surveys, based on the mapping of different 
regions of the sky with a set of contiguous intermediate/narrow 
band filters, provide low resolution photo-spectra for each re¬ 
gion of the sky (hence performing like a low resolution IFU 
with PSF-limited spatial resolution), with the survey depth as the 
only selection criterion and without the typical spectroscopic un¬ 
certainties in the flux calibration. This opens an unprecedented 
way to progress in our understanding about galaxy evolution 
through the study of millions of homogeneous galaxy SEDs, 
both spatially-resolved in the more nearby Universe and inte¬ 
grated. 

This paper is devoted to presenting MUFFIT (MUlti-Filter 
FITting in photometric surveys) a generic code specifically de¬ 
signed for analysing the stellar content of galaxies with available 
multi-filter data (dealing with the technical peculiarities and the 
big amount of high-quality photometric data available in multi¬ 
filter surveys), as well as to show its functionalities, set the accu¬ 
racy and typical uncertainties in the retrieved stellar population 
parameters, and ultimately test it with real data. In this sense, 
we make use of the ALHAMBRA database as a test-bench for 
MUFFIT, not with the aim of performing a thorough stellar pop¬ 
ulation analysis of the galaxies in ALHAMBRA (which consti¬ 
tutes the matter of the next papers of this series), but to compare 
the stellar population results derived from MUFFIT with similar 
studies in the literature, allowing us to assess on its reliability 
and on the feasibility of this kind of techniques to accurately ex¬ 
plore the stellar content of galaxies. 

In the following items the main conclusions of this work are 
summarized; 

- Using as input a set of SSP models that explores different 
stellar population parameters, MUFFIT builds photometric 
predictions of bands at different redshifts and extinctions. 
Eor the present work, the stellar population parameters that 
are considered are just age and metallicity, although in a gen¬ 
eral case even the IME slope and the a-enhancement can also 
be retrieved if SSP models properly account for them. In ad¬ 
dition, the survey photometry is corrected of MW dust ef¬ 
fects, as the colour terms introduced by MW dust may play 
an important role not only in the stellar masses derived us¬ 
ing SED-fitting techniques, but also in the retrieved stellar 
population parameters. 

- MUFFIT compares the multi-filter fluxes of a given galaxy 
with the photometric predictions of a reasonable mixture of 
two SSPs, respectively one younger and one older than the 
mean age provided by a single SSP fitting. The mixture of 
two SSPs determined by the last prior, specific for each in¬ 
dividual galaxy, is a relevant improvement with respect to 
the fitting of one SSP only, since it better represents a pop¬ 
ulation mainly composed by an old component that suffer a 
later burst or star-formation episode. The stellar population 
parameters (in this work the age and metallicity weighted 
by both luminosity and mass, extinction, redshift, and stellar 
mass) provided by MUFFIT are constrained by the use of 
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Fig. 18. As Fig.[T^but using the MIUSCAT SSP models instead of the BC03 ones. 


an error-weighted During the fitting process, MUF¬ 

FIT removes those bands that are affected by emission lines, 
improving the quality of the fitting and restricting the plau¬ 
sible redshift space, as in a general case the redshift of the 
galaxy is treated as another free parameter to be determined. 
MUFFIT is not limited to providing the parameters of the 
best fitting model, but also explores the parameter space us¬ 
ing the proper photometric uncertainties in each band by a 
Monte Carlo method, reinforcing the parameter predictions 
as it provides their statistical uncertainties too. In addition, 
MUFFIT also computes and provides the A:-corrections of 
each galaxy from the same mixture of models at rest-frame. 

- Specifically for the ALHAMBRA data, we study the intrin¬ 
sic uncertainties in redshift, extinction, age, metallicity and 
stellar mass that appear when diagnosed by MUFFIT. Us¬ 
ing the typical distribution of errors for the RS galaxies in 
ALHAMBRA, we construct mock galaxies with an average 
S IN per filter of 20, obtaining typical uncertainties (RMS) of 
CTj ~ 0.01, (TAv ~ 0.11, CTAge ~ 0.10 dex, Cr[Fe/H] ~ 0.16 dex, 
and ~ 0.08 dex. In no case there are systematic errors 
that are statistically significant. 

- Despite MUFFIT is not a generic photo-z code, using as in¬ 
put the redshift PDFs provided by external photo-z codes, 
MUFFIT returns fine-tuned redshift values whose accu¬ 
racy is improved by ~ 10 - 20%. We also obtain that the 
photo-z accuracy reached in ALHAMBRA, (Tnmad ^ 0.009, 
presents a negligible impact on the main stellar population 
parameters retrieved using MUFFIT, being more crucial the 
typical uncertainties in the photometry. 

- We have studied with MUFFIT the age-metallicity- 
extinction degeneracy of the ALHAMBRA data at different 
parameter ranges, and having fixed the IMF slope. The age- 
extinction anti correlation is present in all ranges of age and 
metallicity. However, the well known age-metallicity anti 
correlation may turn into a positive correlation for young 
and/or metal poor populations due to the role of the extinc¬ 
tion in reddening the spectral energy distributions of galax¬ 
ies. 

- The stellar mass predictions provided by MUFFIT for a com¬ 
mon sample of RS galaxies in ALHAMBRA are in won¬ 
derful agreement with the stellar masses computed for the 
same galaxies in COSMOS. The dispersion of the compar¬ 
ison, with an RMS of AlogjgM* ~ 0.15 dex, can be fully 
explained by the intrinsic uncertainties of both methods. 

- MUFFIT offers a reliable way to explore emission lines in 
multi-filter surveys. Using a set of emission line galaxies in 
common between SDSS and ALHAMBRA, we demonstrate 
that the residuals provided by MUFFIT for the filters affected 
by emission lines in ALHAMBRA are correlated with the 
strengths of the main emission lines. 


- The age-metallicity loci provided by MUFFIT for a sample 
of z < 0.22 early-type galaxies in ALHAMBRA at different 
stellar mass bins are in very good agreement with the ones 
determined from SDSS data on the basis of spectroscopic 
diagnostics. When we analyse the stellar content of these 
galaxies in ALHAMBRA using their photometric data and 
MUFFIT, our results point out that the more massive early- 
types 10**M*[Mo]) were formed in an earlier epoch than 
their low-mass counterparts (S 10*°M*[Mo]) with a larger 
content in metals, being these differences AAge ~ 3 Gyr and 
A[Fe/H] ~ 0.3 dex. This result agrees with the "downsizing" 
scenario as well. 

- For a subsample of galaxies in common between ALHAM¬ 

BRA and SDSS, a one-to-one comparison between the red- 
shifts, ages, metallicities, and stellar masses derived spec¬ 
trosco pically for the SDSS data (provided by iGallazzi et al.l 
l2005h and those determined from MUFFIT and ALHAM¬ 
BRA reveal good qualitative agreements in all the parameters 
given the uncertainties of both methods, with typical RMS 
for the distribution of differences between both diagnostics 
of (j-SDSS ^ 0.008, ~ 1.6 Gyr, ~ ^.2 dex, 

and ~ 0.19 dex; hence reinforcing the strengths of 

multi-filter galaxy data and optimized analysis techniques, 
like MUFFIT, to conduct reliable stellar population studies. 
Despite the qualitative agreement between ages, in the sense 
that young (old) spectroscopic ages in SDSS are also found 
to be young (old) photometric ages in ALHAMBRA using 
MUFFIT, there exists a systematic difference of ~ 2 Gyr be¬ 
tween the two samples that is explained by the differences 
of using mixtures of SSPs instead of r-models. Despite the 
good agreement between metallicities, it is noticeable that 
the metallicities of ALHAMBRA galaxies within an aperture 
of < 4" present a better agreement with the spectroscopic 
measurements (A[Fe/H] 0.05 dex) than the galaxies with 
larger apertures (A[Fe/H] $ 0.15 dex), pointing towards the 
possibility that aperture differences between SDSS and AL¬ 
HAMBRA and the existence of metallicity gradients drive 
the observed differences (in agreement with G05). There is 
also a good agreement between stellar masses, with a minor 
shift of ~ 0.18 dex that can be explained by the observed 
offset in age. 

To conclude, we demonstrate that MUFFIT is a reliable 
stellar population code for multi-filter galaxy data, which is 
suited and optimized to analyse the stellar content of galaxies 
in ALHAMBRA-like surveys. This opens a new way to explore 
and address stellar population studies of galaxies with multi¬ 
ple photometric bands or colours, as long as the effective spec¬ 
tral resolution is at least the one of ALHAMBRA, allowing us 
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to accurately extract the stellar content of thousands of galax¬ 
ies at higher redshifts, benefited by the large-number statistics 
in comparison with typical spectroscopic datasets at same red- 
shift. With the arrival of the new generation large-scale multi¬ 
filter surveys, like J-PLUS and J-PAS, codes like MUFFIT will 
contribute greatly to shed light in our understanding about the 
formation and evolution of galaxies. 
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